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ABSTRACT 

A  theory  of  the  surface  tension  of  a  auantum  liquid,  pro¬ 
posed  by  Atkins,  is  developed  further  and  applied  to  liquid 
helium  at  0°K.  On  the  basis  of  the  imperfect-gas  model  developed 
by  Gross  and  Pitaevski  the  free  energy  per  unit  area  arising 
from  the  density  non-uniformity  is  calculated.  This  surface 
free  energy  is  then  identified  with  the  !:intrinsic v  surface 
tension  of  liquid  helium. 

The  same  model  of  the  surface  tension  of  a  quantum  fluid 
is  then  used  to  calculate  the  interfacial  surface  tension  in  a 
phase-separated  He3-He4  system  at  T  =  0°K.  The  result  obtained 
is  then  contrasted  with  the  corresponding  result  based  on  a 
''classical"  model  of  the  surface  tension  of  an  ideal  fluid.  A 
corresponding  experimental  study  should  be  able  to  distinguish 
quantitatively  between  the  two  predictions . 

We  also  discuss  the  effect  of  density  non-uniformities  on 
the  lambda  transition  in  He4  films.  Semi-quantitative  evidence 
is  obtained  that  the  lambda  temperature  depression  in  He4  films 
is,  at  least  in  part,  due  to  the  presence  of  excitations  at  the 
free  surface  of  the  film. 

We  also  obtain  an  estimate  of  the  lambda  temperature  de¬ 
pression  in  rotating  helium  as  a  function  of  both  the  angular 
speed  of  rotation  and  the  physical  dimensions  of  the  container. 
The  result  obtained  is  more  general  than,  and  under  special 
conditions  reduces  to,  the  previous  results  obtained  by 
Andronikashvili ,  Mamaladze  and  co-workers. 
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CHAPTER  1 


INTRODUCTION 

The  study  of  liquid  helium  is  of  considerable  physical 
interest  because  it  is  the  only  liquid  that  exhibits  quantum 
effects  on  a  macroscopic  scale.  Before  we  delve  Into  the  study 
of  certain  specialized  problems  concerning  this  liquid,  it 
appears  worthwhile  to  give  a  brief  outline  of  some  of  Its 
important  properties  that  arise  by  virtue  of  these  quantum 
effects . 

Kamerlingh  Onnes  (1)  was  the  first  to  successfully 
accomplish  in  1908  the  liquefaction  of  helium.  He  discovered 
that  the  critical  temperature  of  heliurrr ;was  about  5.2°K,  and 
its  normal  boiling  point  about  4.2°K.  Solid  helium,  however, 
could  not  be  obtained  simply  by  cooling  the  liquid.  Owing  to 
the  large  zero-point  energy  of  the  helium  atoms  (which  arises 
from  their  light  mass)  and  the  relatively  weak  interatomic 
attractive  forces  (which  are  so  because  of  the  closed  structure 
of  the  helium  atoms),  the  solid  phase  can  be  obtained  only  by 
applying  an  external  pressure  of  at  least  25  atmospheres  even 
at  T  =  0°K . 

As  early  as  1911?  Kamerlingh  Onnes  (2)  noticed  that  strange 
things  happened  to  liquid  heliumrwhen  it  is  cooled  below  about 
2.2°K.  As  the  temperature  is  decreased  below  this  value,  the 
liquid,  instead  of  continuing  to  contract,  began  to  expand. 

This  suggested  to  Keesom  and  Wolfke  (3)  the  possibility  of 
helium  undergoing  an  allotropic  modification  at  T  -  2.2°K;  they 
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designated  the  liquid  at  lower  temperatures  as  Hell  and  that  at 
highfer  temperatures  as  Hel . 

The  discovery  that  gave  this  transition  the  name  X-trans- 
ition  was  due  to  Keesom  and  Clusius  (4)  who  found  in  1932 
that  while  the  phase  transition  at  2.l8°K  was  not  accompanied 
by  a  latent  heat,  it  was  indeed  accompanied  by  a  singularity 
in  the  nature  of  the  specific  heat  curve,  whose  shape  resembled 
the  greek  letter  lambda.  Recent  measurements ,’ -by  Fairbank, 
Buckingham  and  Kellers  (5)>  of  the  specific  heat  at  saturated 
pressure  have  shown  that,  as  the  temperature  approaches  T,  c  , 

A  7  S  a,  U 

diverges  in  the  manner  of  the  logarithm  of  the  magnitude  of 
the  temperature  interval  |  AT^|  =  |T-T^|.  For  the  temperature 

interval  10  <  ] AT |  <  10  ,  the  specific  heat  satisfied  the 

empirical  formula: 


where 

there 


°sat  =  4'55  "  3'00  loSlO 

6=0  for  T  >  T^ ,  and  6  = 
occurs  a  discontinuity  in 
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j  oule 
gm  °K  , 


(1.1) 


1  for  T  <  T^  .  It  is  clear  that 
the  specific  heat  of  the  amount 


A  c 


sat 


=  c 


sat 


(T>TX)  - 


csat(T<V  =  "5-20 


j  oule 

gm  °  K  .  (112) 


A  closely  related  property  of  liquid  helium  was  discovered 
by  Keesom  and  Miss  Keesom  in  1936.  They  noticed  that  as  liquid 
helium  passes  through  the  X-point  into  the  helium  II  phase  there 
occurred  an  enormous  increase  in  its  thermal  conductivity . 

Earlier  in  1930,  Keesom  and  van  den  Ende  (6)  discovered, 
rather  accidentally,  one  of  the  most  important  properties  of 
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liquid  helium.  They  observed  that  liquid  helium,  when  cooled 
below  2.l8°K,  could  manage  to  escape  through  extremely  tiny 
leaks  in  the  apparatus  which  had  been  completely  leakproof  at 
higher  temperatures  to  both  helium  I  and  gaseous  helium.  This 
observation  seemed  to  indicate  an  enormous  drop  in  the  vis¬ 
cosity  of  the  liquid  on  passing  down  through  the  lambda  point. 

Early  measurements  on  the  viscosity  of  liquid  helium, 
however,  yielded  conflicting  results.  Experiments  based  on  the 
oscillating  disc  method  (7)  showed  that  the  viscosity  of 
helium  II,  though  it  decreases  continually  with  decreasing 
temperature,  is  not  very  different  from  the  viscosity  of 
helium  I.  On  the  other  hand,  experiments  done  by  Kapitza  (8) 
and  by  Allen  and  Misener  (9)  using  the  capillary  flow  method, 
seemed  to  indicate  that  the  viscosity  of  helium  II  was  several 
orders  of  magnitude  lower  than  that  of  helium  I. 

Other  properties  which  indicate  the  inadequacy  of  a 
classical  description  for  liquid  helium  were  the  so-called 
fountain  effect  and  the  existence  of  supersurface  films .  Both 
these  properties  seemed  to  indicate  that  the  heat  transfer  and 
the  mass  transfer  in  helium  II  were  inseparably  interconnected  (10). 

The  first  to  tackle  the  problem  theoretically  was 
F.  London  (11)  who  in  1938  advanced  the  hypothesis  that  the 
phenomenon  of  Bose-Einstein  condensation,  with  due  allowance 
made  for  the  presence  of  interparticle  interactions  and  for  the 
liquid  state  of  the  substance,  might  be  the  essential  cause  of 
the  X-transition  in  liquid  helium.  This  is  made  highly  plausible 
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by  the  simple  fact  that  whereas  the  heavier  Isotope  of  helium 
(namely  He4)  displays  the  phenomenon  of  A-transition  the  lighter 
isotope  (namely  He3)  does  not.  Moreover,  the  temperature  at 
which  the  ideal  Bose  gas,  with  the  same  atomic  mass  and  density 
as  liquid  He4,  undergoes  Bose-Einstein  condensation  turns 
out  to  be  about  3-l4°K,  which  is  of  the  very  order  of  magni¬ 
tude  as  the  transition  temperature  .in  liquid  helium.  London 
went  on  to  suggest  that  the  superflow  of  Hell  might  be  supposed 
to  be  connected  especially  with  the  condensed  phase  of  the  Bose 
gas.  Further,  there  is  some  similarity  (though  not  very  much) 
in  the  temperature  dependence  of  the  specific  heat  in  the  two 
cases.  Whereas  liquid  helium  shows  a  logarithmic  discontinuity 
in  the  specific  heat  at  the  transition  point,  the  ideal  Bose 

gas  shows  a  discontinuity  only  in  the  derivative  of  the 

\9T  /V 

specific  heat.  Of  course,  in  the  Bose  gas  model,  the  specific 

heat  C  does  become  infinite  at  the  transition  temperature. 

]? 

An  important  difference  is  met  with  in  connection  with  the 
behaviour  of  at  l°w  temperatures.  While  liquid  helium 

shows  a  T3-behaviour ,  the  ideal  Bose  gas  shows  a  T3/2  behaviour. 

Despite  these  shortcomings,  London’s  hypothesis  did 
provide  a  theoretical  basis  for  the  two-fluid  model  of  liquid 
helium,  which  had  been  suggested  by  Tisza  (12)  on  purely  em¬ 
pirical  grounds.  This  model  rested  on  the  following  major 
assumptions  (13): 
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1)  There  exists  in  liquid  helium  II  a  superfluid  com¬ 
ponent  plus  an  excited  normal  fluid  component.  At  T  =  0°K, 
the  whole  of  the  liquid  is  superfluid.  As  the  temperature  is 
raised  excitations  are  created  whose  density  increases  mono- 
tonically  with  temperature;  these  excitations  constitute  the 
normal  fluid.  At  T^T^ ,  the  whole  of  the  liquid  is  in  the  normal 
state . 

2)  The  viscosity  and  the  entropy  of  the  superfluid  com¬ 
ponent  are  identically  zero.  The  viscosity  and  the  entropy  of 
liquid  helium  II  are,  therefore,  due  to  the  normal  component 
alone . 

These  assumptions  at  once  remove  the  apparent  discrepancy 
in  the  viscosity  measurements.  In  the  capillary  flow  experi¬ 
ments  the  normal  component  is  completely  held  back  and  the  flow 
is  all  due  to  the  superfluid  component.  Consequently  the 
measured  viscosity  is  practically  zero.  In  the  oscillating 
disc  experiments,  however,  the  superfluid  component  remains 
unattached  while  the  normal  component  becomes  Involved  in  the 
oscillatory  motion.  The  measured  coefficient  of  viscosity  then 
depends  on  the  fraction  PR/p  of  the  normal  component  of  the 
liquid . 

The  explanation  for  the  existence  of  second  sound  In 
liquid  helium  also  follows  easily  from  these  assumptions.  If 
the  temperature  at  some  point  in  the  fluid  is  suddenly  raised, 
excitations  will  be  created  at  that  point  leading  to  a  non¬ 
uniformity  in  the  density  of  the  normal  component.  A  resulting 
flow  of  normal  fluid  away  from  that  point  will  occur  with  a 
corresponding  flow  of  superfluid  component  towards  it.  Entropy 
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non-uniformities  can  thus  be  propagated  as  temperature  waves 
through  liquid  helium  II.  This  phenomenon  was  first  studied 
by  Tisza  (12.)  in  1940  and  Landau  (14)  in  1941. 

Although  experiments  bore  out  many  of  the  predictions  of 
the  two-fluid  model,  it  remains  primarily  an  intuitive  model, 
not  based  on  any  rigorous  theory  except  on  the  relevance  of 
the  Bose-Einstein  condensation  to  the  observed  phenomena  of 


superfluidity . 

The  two-fluid  model 

was  given  a  new  life  by  Landau  ( 1 4 ) . 

In  1941,  Landau  suggested 

that  the  liquid  helium  be  regarded  as 

a  quasi-continuum,  whose  low-lying  excitations  were  quantized 
in  the  same  manner  as  Debye  modes  in  a  solid.  Landau  assumed 
two  kinds  of  elementary  excitations:  (l)  phonons,  which  are 
the  quanta  of  longitudinal  sound  waves,  governed  by  the  energy 
spectrum  £=pc,  where  p  is  the  momentum  of  the  phonon  and  c  is 
the  velocity  of  sound  in  the  liquid,  and  (ii)  rotons,  whose 
energy  is  assumed  to  be  separated  from  the  lowest  state  by  an 
energy  gap  A.  Both  types  of  excitations  are  supposed  to  have 
an  effective  mass  and  momentum  which  are  ultimately  attributed 
to  the  normal  component  of  the  liquid.  The  energy  spectrum  of 


the  rotons  may  be  written 

as 

E  =  A 

.  <P  -  Po>2 

2y  ,  (1.3) 

where  A,  pQ  and  y  are  the 

parameters  of  the  spectrum,  The 

precise  values  of  these  parameters  are  hard  to  obtain  but  have 
been  found  by  means  of  neutron  scattering  experiments  (1-5);  thus 
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we  have 

A/k  =  8.65°K,  P0/ta  =  1.91A  \  y  =  0.l6m^e. 

1  The  existence  of  another  kind  of  excitation  has  been 
proposed  by  Kuper  (16);  these  are  associated  with  the  vibrational 
modes  of  the  free  surface  of  the  liquid  and  are  considered  to 
be  of  the  nature  of  "quantized  capillary  waves"  or  the  so-called 
ripplons .  The  energy  spectrum  of  these  excitations  is 

go2  =  g  k  +  ,  (1.4) 

where  the  first  term  arises  from  the  influence  of  gravity  on 
the  vibrations  of  the  surface  while  the  second  term  arises 
from  the  surface  tension  *of  the  liquid.  The  concept  of  quantized 
surface  modes  of  vibration  has  proved  very  useful  in  under¬ 
standing  the  temperature  variation  of  the  surface  tension  of 
liquid  helium  and  also  leads  to  a  significant  contribution 
towards  the  zero-point  surface  tension  of  the  liquid.  Of  course, 
the  property  of  surface  tension  arises  intrinsically  from  the 
interparticle  interactions  in  the  system.  A  study  along  these 
lines  has  been  carried  out  in  Chapters  2  and  3 • 

Other  studies  based  on  the  two-fluid  model  of  helium  II 
that  have  been  carried  out  here  relate  to  the  lambda  temperature 
depression  that  occurs  in  helium  films  and  is  supposed  to  occur 
in  rotating  helium.  These  problems  have  been  discussed  in 
detail  in  Chapters  5  and  6 . 
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A  number  of  microscopic  theories  have  also  been;  .developed 
in  order  to  understand  the  physical  behaviour  of  liquid  helium 
from  a  more  fundamental  standpoint .  Difficulties  in  developing 
a  workable  theory  from  first  principles,  however,  stem  from  the 
fact  that  the  substance  under  study  is  in  the  liquid  state, 
which  of  all  the  states  of  matter  is  the  least  understood. 
Nevertheless,  the  theories  of  Bogolyubov( 17)  and  Lee,  Huang, 
and  Yang  (18)  have  succeeded  in  explaining  some  of  the  quali¬ 
tative  features  of  the  liquid  helium  behaviour.  These  theories 
have,  almost  invariably,  employed  a  hard-sphere  interaction 
as  a  model  potential  between  the  helium  atoms.  However,  in 
order  to  study  liquid  helium  on  a  more  realistic  basis, 
especially  in  regard  to  properties  that  belong  exclusively  to 
the, liquid  state,  one  must  not  overlook  the  attractive  part  of 
the  interaction.  Without  this  part  of  the  interaction,  the 
system  cannot  even  exist  in  the  liquid  state.  The  attraction 
may,  of  course,  be  too  weak  to  cause  two-body  bound  states  in 
the  medium,  but  it  must  at  any  rate  be  strong  enough  to  form 
an  N-body  bound  state.  The  N-body  system  will  consequently 
possess  an  equilibrium  density  at  which  the  total  energy  of 
the  system  is  negative.  'Therefore,  when  attractive  inter¬ 
actions  of  reasonable  strength  are  present  it  is  no  longer 
necessary  to  fix  the  volume  of  the  system  by  means  of  external 
boundaries.  In  fact,  if  the  volume  of  the  container  is  larger 

than  that  determined  by  the  equilibrium  density  of  the  liquid, 
./Surface 

a  free /appears .  A  study  of  the  nature  of  this  free  surface 
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and  a  calculation  of  the  free  energy  associated  with  it,  on  the 
basis  of  the  imperfect-gas  model  of  Gross  and  Pitaevski, 
constitutes  the  first  major  content  of  this  thesis.  By  means 
of  this  approach,  we  obtain  a  theoretical  estimate  of  the 
"intrinsic”  surface  tension  of  liquid  helium;  see  Chapter  3. 

We  propose  to  review  other,  more  particular,  properties 
and  concepts  relating  to  liquid  helium  in  the  introductory 
sections  of  the  relevant  chapters  of  the  thesis. 
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CHAPTER  2 


TEMPERATURE  DEPENDENCE  OF  THE  SURFACE  TENSION  OF  HELIUM  II 
2.1  INTRODUCTION 

The  surface  tension  o.f  liquid  helium  was  originally  measured 
by  van  Urk,  Keesom,  and  Kamerlingh  Onnes  (19)  in  1925  and  by 
Allen  and  Misener  (20)  in  1938.  Extrapolating  the  results 

obtained  a,t  0°K,  they  observed  that,  starting  from  a  value  of 

_2 

0.35  erg  cm  ,  the  surface  tension  at  first  decreases  gradually 

but  begins  to  decrease  more  ,  rapidly  as  the  critical  point  is 

approached.  Recent  experimental  work  by  Atkins  and  ltfarahara  (21) 

4 

on  the  determination  of  the  surface  tension  of  helium  for 
temperatures  down  to  0.35°K  has  shown  that  the  temperature  depen¬ 
dence  of  the  surface  tension  can  be  expressed  by  an  empirical 
relationship  of  the  form 

a(T)  =  a ( 0 )  -  aTn  ;  (2.1) 

here,  a(0)  is  the  limiting  value  of  the  surface  tension  as 
T-*-OdK,  while  a  and  n  are  adjustable  parameters.  The  best  fit, 
jising  the  method  of  least  squares,  was  obtained  with  a(0)  = 

0.3729  e.rg  cm-2,  a  =  0.0081  erg  cm-2  0K_n,  and  n  =  2.5±0.2.  A 
complete  theoretical  understanding  of  the  observed  behaviour  of 
a(T),  as  depicted  by  Eq .  (2.1)  was,  however,  still  lacking. 

The  most  satisfactory  theory  we  have  so  far  of  the 
surface  tension  of  liquid  helium  II  is  the  one  due  to  Atkins 
who,  in  an  earlier  paper  (22),  suggested  that  the  major  con¬ 
tribution  to  the  surface  energy  of  this  liquid  comes  from  the 
quantized  surface  modes  of  vibration  known  as  ripplons .  Atkins 
assumed  these  surface  modes  to  be  similar  in  nature  to  the 
macroscopic  capillary  waves  which  have  a  frequency-dependent 
phase  velocity  given  by  (23) 
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2.2 


( 27rav/p ) 


1  /  3 
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(2.2) 


where  a  is  the  surface  tension  of  the  liquid,  p  is  its  density 
and  v  is  the  frequency  of  the  wave.  In  the  spirit  of  the 
Debye  theory  of  solids,  Atkins  chose  a  cutoff  frequency  v  , 
such  that  the  total  number  of  normal  modes  in  the  surface  was 
equal  to  the  total  number  of  atoms  in  a  monomolecular  layer 
at  the  surface.  This  fixed  the  cutoff  frequency  at  about  1011 
per  second  and  the  corresponding  characteristic  temperature 
0  ( =h\)  /k)  at  a  few  degrees  Kelvin. 

G  C 

Thus,  according  to  Atkins’  theory,  a  part  of  the  observed 
surface  tension  a(T)  of  liquid  helium  arises  from  the  presence 
of  quantized  surface  modes,  which  themselves  arise  due  to  the 
surface  tension  of  the  liquid.  Accordingly,  we  must  invoke 
an  ’’intrinsic”  surface  tension,  ck(T)  say,  of  this  liquid, 
with  the  result  that  the  observed  surface  tension  cr(T)  can  be 
expressed  as 


(2.3) 


a(T)  =  a.(T)  +  a 

1  V' 


being  the  contribution  arising  from  the  vibrational  modes 
of  the  surface  which  would  depend  upon  the  temperature  of  the 

liquid  explicitly  as  well  as  implicitly - the  latter  through 

the  temperature  dependence  of  a^(T)  itself.  Of  course,  it  will 
be  natural  for  us  to  look  for  the  origin  of  ck(T)  in  the  basic 
interactions  among  the  atoms  of  the  fluid.  This  problem  will 
be  pursued  in  Chapter  3- 


In  the  limit  of  0°K,  Eq .  (2  .3)  reduces  to  (22) 


(2.4) 
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2.3 


the  second  term  on  the  right  representing  the  zero-point  con¬ 
tribution  of  the  surface  modes .  In  order  to  appreciate  the 
relative  magnitudes  of  the  two  terms  in  (2.4),  we  note  that  to 

obtain  the  experimental  estimate  of  0.37  erg  cm  for  a£0),  we 

_  o 

require  ck(0)  to  be  about  0.14  erg  cm  ,  leaving  about  0.23 

_2 

erg  cm  for  a  (0).  The  contribution  of  the  surface  modes  is, 
therefore,  more  than  appreciable! 

At  non-zero  temperatures,  especially  for  T<<0c,  the 
added  contribution  of  the  surface  modes  (as  arising  from  the 
explicit  dependence  on  temperature)  turns  out  to  be  (4) 

_1'55  (jCoxT  3^75-  =  ~6  •  5  x  1CT3  T7/3  erg  cm-2  .  (2.5) 

The  exponent  7/3  is  quite  consistent  with  the  empirical  value 
of  2.5  ±  0.2’  however,  the  numerical  factor  here  Is  about  20$ 
lower  than  the  corresponding  empirical  value.  It  Is,  however, 
quite  imaginable  that  the  temperature  dependence  of  a^(T) 
itself  could,  to  some  extent,  account  for  the  remainder. 

An  alternative  explanation  of  the  temperature  dependence 
of  the  surface  tension  of  liquid  helium  has  been  attempted  by 
Singh  (24),  who  adopted  the  Ideal  gas  model  for  this  liquid 
but  made  use  of  a  better  enumeration  of  the  eigenfunctions  in 
a  bounded  continuum,  viz. 

g(p)  dp  =  P2dp  +  jjW  pdp  +  .  .  .  ,  (2.6) 

where  V  is  the  volume  of  the  enclosure,  and  A  Its  surface  area. 
The  positive  sign  of  the  surface  term  is  due  to  the  fact  that 
Neumann  boundary  conditions  have  been  adopted,  i.e.  ^  =  0  at 
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the  surface.  The  above  enumeration  takes  into  account  the 


so-called  surface  effects  in  the  expression  for  the  density  of 
states .  Singh  thereby  obtained  for  the  temperature-dependent 
part  of  the  surface  tension  the  expression 


— tt  m  -k2  C  ( 2 )  T: 

2k7 


•7.5x10  erg  cm  ^ 


(2.7) 


here,  m  is  the  mass  of  a  helium  atom,  while  the  Riemann  function 
C(2)  is  equal  to  tt .  Comparing  (2.7)  with  the  empirical  result 
of  Atkins  and  Narahara,  one  finds  that  the  numerical  factor  is 
now  much  closer  to  the  corresponding  empirical  value,  but  that 
the  exponent  is  significantly  low.  Prom  the  point  of  view  of 
agreement  with  the  experimental  data,  therefore,  there  is  not 
much  reason  to  choose  one  or  the  other  of  the  two  alternatives 
(2.5)  and  (2.7).  Prom  a  physical  point  of  view,  however,  we 
note  that  whereas  the  theoretical  basis  for  (2.5)  is  suffi¬ 
ciently  sound,  the  adequacy  of  using  the  ideal  gas  model  for 
liquid  helium,  which  led  to  (2.7)  is  hardly  acceptable.  Of 
course,  Singh’s  argument  does  advocate  a  closer  examination  of 
the  surface  effects  than  has  hitherto  been  done;  however,  in 
studying  the  contribution  of  these  effects  to  the  density  of 
states  and  hence  to  the  free  energy  of  liquid  helium,  we  must 
consider  a  gas  of  actual  excitations,  namely  phonons  and  rotons, 
rather  than  one  of  non-interacting  atoms .  In  the  absence  of 
interactions  there  would  be  no  basis  for  expecting  0^(0),  and 
hence  a(0)  to  be  nonzero  at  T=0°K. 
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2.5 


In  this  chapter  we  shall  settle  the  issue  in  regard  to 
the  surface  effects  that  arise  from  a  better  enumeration  of 
the  density  of  states  in  a  gas  of  excitations.  In  Sec.  2.2, 
we  shall  consider  the  thermodynamics  of  an  n-dimensional  ideal 
Bose  gas  of  quasi-particles  (excitations)  taking  into  account 
"surface"  effects  that  may  be  important.  In  Sec.  2.3  the 
surface  energy  contributions  made  by  phonons  and  rotons  are 
calculated  and  it  will  be  shown  t'hat ,  regardless  of  boundary 
conditions  chosen,  these  are  negligible  compared  to  the  free 
energy  contributions  due  to  ripplons.  The  temperature 
dependence  of  the  surface  tension  of  helium  II  can  thus  be 
adequately  explained  by  the  existence  of  quantized  sur-pace 
modes  of  vibration. 


. 
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2.2  THE  n -DIMENSIONAL  BOSE  GAS 

In  order  to  study  /the  physical  properties  of  a  statistical 
system  a  knowledge  of  the  asymptotic  distribution  of  single 
particle  states,  which  are  eigenfunctions  of  the  Schrodinger 
equation 


(2.8) 


is  essential.  Customarily,  the  number  of  eigenstates  g(p) 
with  momenta  not  exceeding  a  particular  value  p  Is  approxi¬ 
mated  by  the  familiar  Rayleigh-Weyl  expression 


g(p)  =  Pp  (p/j?)3 


(2.9) 


where  V  is  the  volume  of  the  enclosure.  Of  course,  the 
validity  of  this  expression  rests  on  the  assumption  that  the 
linear  dimensions  of  the  enclosure  are  much  larger  than  the 
characteristic  wavelength  (,h/p)  .  A  more  rigorous  study  of  the 
density  of-  states,  as  allowed  by  (2.8)  shows  that  the  second 
term  of  the  asymptotic  expression  for  g(p),  of  which  (2.9)  is 
the  main  term  is  directly  proportional  to  the  surface  area  S 
of  the  enclosure  (25)  (26).  This  surface  term  is  naturally  of 
a  very  direct  physical  interest  to  us  for,  when  taken  into 
account  in  calculating  the  extensive  properties  of  a  physical 
system,  it  leads  to  what  might  be  called  a  "surface  contri¬ 
bution"  to  the  customary  bulk  values  of  these  properties . 


In  order  to  appreciate  the  accuracy  made  possible  by  the 
inclusion  of  the  surface  term,  Pathria  (25)  has  tabulated  a  set 
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numbers  adapted  from  an  earlier  work  by  Bolt  (27) }  who  computed 
the  actual  frequency  distribution  of  acoustical  eigentones 
(under  the  Neumann  boundary  condition)  for  a  number  of  rooms 
of  various  sizes.  In  a  typical  case,  when  the  room  size  was 
10’  by  15’  by  30’,  the  numerical  value  of  g(v), taking 
c-1125  ft/sec,  were  as  given  in  the  last  column  of  Table  2.1. 

Table  2.1 


Frequency 

V 

(cps) 

g(v)  asymptotic 

g(v)  actual 

(from  Bolt) 

Volume  term 

?jlone 

Volume  term 
and 

Surface  term 

60 

3 

7 

8 

100 

13 

24 

25 

150 

^5 

70 

74 

200 

106 

151 

159 

250 

207 

277 

283 

300 

358 

458 

464 

The  second  column  of  the  table  gives  the  numbers  due  to  the 
volume  term  alone,  while  the  third  column  takes  into  account  the 
surface  term  as  well.  Clearly,  the  surface  term  makes  a  very 
important  contribution  to  g  (v)  ;  at  the  same  time,  the  final 
remainder  appears  to  be  very  small. 

For  the  sake  of  complete  generality,  the  enclosure  that  we 
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shall  consider  is  an  n-dimensional  domain,  say.  For  such 
a  domain  Fedoso.v  (26)  has  established  the  following  asymptotic 
formula : 


jf(K-k)C  dg(k) 
o 


cm 
TT+i)  ! 


K^+1  +  0(Kn  X), 


(2.10) 


so  that,  with  £  =  0,  one  obtains 

g ( K )  =  i  a,  K1  +  0 ( Kn_1 ) . 
1 


(2.11) 


The  coefficients  a^  here  are  related  to  certain  geometrical 
measures  of  the  domain,  viz. 


(2.12) 


and 


L 

n-1 


£  mes  ,  -.x  S  . 

•  (n-1)  j 

+  _n+l 


(n-1) 


(A 


(2.13) 


where  Sj ^  is  a  t-dimensional  "face"  of  the  domain;  the  upper 
sign  corresponds  to  the  Dirichlet  boundary  condition  ip  =  0 , 
and  the  lower  sign  corresponds  to  the  Neumann  boundary 

condition  =  0. 
dn 


For  a  three-dimensional  domain  of  arbitrary  shape,  Eq .  (2.12) 
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and  (2.13)  yield  the  formula 


(2. 11’ ) 


Now,  the  error  term  as  such  rather  appears  to  vitiate  the  role 
of  the  surface  term  in  the  expression  for  g(K) .  However, 
Pathria  (25)  has  provided  sufficient  computational  evidence 
(see  Table  2.1)  to  show  that  the  error  term  in  (2.11f)  is  in 
effect  much  inferior  to  the  surface  term,  so  that  Eq .  (2.11’) 
may  be  replaced  by 


(2.11?!') 


Using  (2.12)  and  (2.13)  we  can  calculate  the  thermo¬ 
dynamic  functions  of  an  n-dimensional  Bose  gas.  We  have,  quite 
generally , 


=  -  Z  In  (1  -  z  e  3e£-)  . 

£ 


(2.14) 


The  summation  can  be  converted  into  an  integral  employing  the 
expression  (2.11),  along  with  (2.12)  and  (2.13),  for  the 
density  of  states.  Again,  for  the  sake  of  generality,  we  take 
the  energy-momentum  relationship  for  the  quasi-particles  to  be 


E  =  aps .  We  then  obtain 
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PV  =  k  T 


In  2.= 


7T 


') 


(z)  V, 
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(kT) 


+  l 


n 
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"A 


(2.15) 


where  g  (z)  =  ? 

n  £=1 


The  internal  energy  can  be  found  from  Eq .  (2.14)  yielding 
the  result 


The  number  of  particles  in  the  excited  states  can  also  be 


found  with  the  result 
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TT 


n/2  fn 


(-\  !  g  (z)  V 
\s,/  6n  v  '  n 


hn  an/s 


CkT) 


n/s 


(2.17) 


From  the  last  relation  it  is  easy  to  show  in  general  the 
conditions  necessary  For  Bose-Einstein  condensation  to  occur 
in  an  infinite  n-dimensional  Bose  gas.  In  face,  true  Bose- 
Einstein  condensation  can  occur  only  if  the  expression  for  the 
total  number  of  excited  particles  (2.17)  is  bounded.  Then 
if  the  total  number  of  particles  is  greater  that  the  maximum 
number  of  particles  the  excited  states  can  hold,  i.e.  If 
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then  the  ground  state  of  the  system  will  be  occupied.  The 

transition  temperature  T  is  then  defined  as 

c 


(2.18) 
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A  non-zero  transition  temperature  is  possible  if  and  only  if 

the  Riemann  zeta  function  c(  — )  remains  bounded.  This  is  true 

s 

only  if  n>s.  Hence  for  an  infinite  three  dimensional  Bose  gas, 
since  n  =  3  and  s  =  2,  a  true  Bose-Einstein  condensation  temp¬ 
erature  does  exist.  In  the  case  of  an  infinite  two-dimensional 
Bose  gas  Bose-Einstein  condensation  does  not  occur  since  C(l) 
is  unbounded  and  hence  all  the  particles  can  be  accommodated 
in  the  excited  states.  Peculiarly  enough,  for  an  extreme 
relativistic,  infinite  two-dimensional  Bose  gas,  since  a  =  pc 
(where  c  Is  the  speed  of  light)  true  Bose-Einstein  condensation 
does  occur  since  now  n  =  2,  s  =  1  and  c(2)  is  bounded. 

Now,  since  F  =  yN  -  PV,  where  y  is  the  chemical  potential 
of  the  Bose  gas,  then  below  the  transition  temperature  (where 
y  =  0)  we  have  P  =  -PV.  Considering  only  the  surface  term  of 
(2.15),  we  have,  in  the  case  of  ideal-gas  particles  (n  =  3, 
a  =  •|m,  s  =  2)  , 

Ps/A  =  ±  (kT)2  =  ±7-5  X  10  “3T  erg  cm-2  >  (2.19) 

which  Is  exactly  the  result  obtained  by  Singh.  Eq .  (2.19)  can 
also  be  derived  without  making  use  of  the  explicit  expression 
for  g(K).  In  that  case,  one  works  with  the  summations  as  such 
by  applying  the  theta-function  transformation  to  the  sums;  for 
the  details  see  Appendix  I.  Still  another  way  of  deriving  this 
result  is  through  the  construction  of  the  Mellin  transforms  of 
the  relevant  summands  and  evaluating  the  resulting  complex 
contour  integrals  (28).  For  the  phonons  (n  =  3,  a  =  c  (the 
speed  of  sound),  s  =  1)  one  obtains 
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2.13 


F  /A 

5! 


+ 


TT  C(3) 

2h  V2 


(kT)3  =  ±  1.96 


_4  3 

x  10 


erg  cm 


(2.20) 


This  result  can  also  be  derived  by  the  use  of  Mellin  trans¬ 
forms  ( 28 )  . 

For  rotons  it  is  necessary  to  make  a  separate  calculation 

because  their  energy-momentum  relationship  does  not  fall  into 

/  \  2 

the  category  E  =  aps .  One  has,  in  this  case,  E  =  A+  ^p~Pgr<--  . 

A  straight-forward  calculation  making  use  of  formulae  (2.11?) 
gives 


co 


F  /A 
s 


7T 


1  kT  ’2hr  J  6 
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A+(P-PC)2^ 


2y 

-A/kT 


/kT 
P  dp 


2h' 


Po  e 


( kT )  3/2 


_p  -A/kT  _p 

±3.1x10  e  T3/2  erg  cm 
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2.3  SURFACE  FREE  ENERGY  OF  RIPPLONS 

Before  applying  the  thermodynamics  of  Section  (2.2)  to 
calculate  the  surface  free  energy  of  the  ripplons  we  shall  have 
to  derive  the  necessary  energy-momentum  relation  for  surface 
waves . 

The  basic  equations  required  in  order  to  derive  the  dis¬ 
persion  relation  for  hydrodynamic  capillary  waves  are  (i) 
Laplace’s  equation  (29) 


(2.22) 


where  C  is  the  normal  displacement  of  the  surface,  and  (ii) 
Euler's  equation  for  an  Incompressible,  irrotational  fluid  in 
a  gravitational  field 


^  grad  v2  =  -  grad  (|-  +  gz) 


(2.23) 


In  this  case  we  have  potential  flow  and  can  write  v  =  grad  <J>, 
so  that  (j)  is  the  velocity  potential.  Equation  (2.23)  then 
reads 


grad  (||-  +  v2  +  P/p  +  gz)  =  0 


and  hence,  choosing  the  arbitrary  function  _p(t)  to  be  zero. 


P 


(2.24) 


Considering  only  small  velocities  of  the  fluid,  and 
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neglecting  the  term  75-  pv  ,  we  substitute  Eq .  (2.21)  Into 
Eq .  (2.22)  and  obtain  the  relation 


,  3  d) T  ,  3  2  r  ,  3 2  C x 

pgz  +  - o  (t^tt  +  T^rr)  =  0  . 


3x' 


3y 


(2.25) 


Here  we  have  redefined  the  potential  6  T  =  <J>  +  (;PQ/p)t;  this 
makes  no  essential  difference  since  v  =  grad  (p  =  grad  6  f  • 

We  can  differentiate  this  relation  with  respect  to  t  and 
replace  ^  by  obtaining 


34>  ,  924>  _  9  /924)  .  32<K 

pg  3z  +  p  Ft*  a  3  z  +  3y^  0> 


(2.26) 


Assuming  now  a  plane  wave  propagating  in  the  x-direction, 
say  (p  =  A  e  cos  (kx-wt)  and  substituting  into  (2.26)  we 
obtain  the  dispersion  relation 


w2  =  gk  + 


crk: 


(2.27) 


For  long  wavelengths,  such  that  k<</gp/a ,  we  have  pure  gravity 
waves .  In  the  opposite  case  of  short  wavelengths  we  have 
capillary  waves,  or  ripplons 


w‘ 


ak3 

P 


(2.28) 


We  can  now  use  Equation  (2'..15)  to  calculate  the  free 
energy  of t bese  excitations  and  consequently  the  contribution  to 
the  surface  tension  of  helium  II  below  the  transition  point. 
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For  ripplons,  the  energy-momentum  relationship,  when  the 
effect  of  the  gravitational  field  can  be  neglected,  is 


E  = 


(— ) 

vph; 


V2 


3/2 

P 


(2.29) 


O’  7 

so  that  we  have  a  =  (— r)  ,  s 

stitute  a  surface  phenomenon 
some  bulk  value),  the  desired 
from  the  main  term  in  (2.15) 
thus  obtains 


=  3/2. 
( rather 
result 
with  n 


Since  ripplons  con- 
than  a  correction  to 
in  this  case  follows 
2  and  Vp  =  A .  One 


2/3 


7/3 


p/a  =  C&j)  h  (^)  r(j)  Z(l) 


=  -6.5  x  10  ^  T7/3  erg  cm 


-2 


(2.30) 


This  is  precisely  the  result  obtained  earlier  by  Atkins  (22). 

A  comparison  of  (2.20)  with  (2.30)  readily  shows  that 
the  surface  effect  of  the  phonon  component  is  hardly  a  few 
percent  of  the  direct  contribution  of  the  ripplons .  The 
surface  effect  of  the  roton  component,  as  given  by  (2.21),  is 
still  less  by  a  factor  of  10  or  more.  Moreover,  a  bh.sic 
difficulty  with  the  surface  effects  (2.20)  and  (2.21)  is  that 
they  are  so  strongly  dependent  on  the  nature  of  the  boundary 
conditions  used.  Thus,  we  see.  that  the  observed  temperature 
dependence  of  the  surface  tension  of  liquid  helium  is,  in  the 
main,  given  by  the  direct  surface  energy  of  the  quantized  waves 
on  the  surface  of  the  liquid. 
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CHAPTER  3 


INTRINSIC  SURFACE  TENSION  OF  HELIUM  II  AT  T=0°K 
3 . 1  INTRODUCTION 

In  this  chapter  we  attempt  to  explain  the  existence  of  the 
limiting  surface  tension  a(o)  on  the  basis  of  the  atomic  inter¬ 
actions.  The  theoretical  model  employed  for  this  investigation 
consists  of  a  dilute,  imperfect  gas  of  bosons,  interacting 
through  a  two-body  potential  V(r-r’)  which  can  be  replaced  by 
a  pseudopotential  V^Cr-r’)  where  V*  is  a  constant  that  depends 
upon  the  strength  of  the  interaction.  It  Is  wall  known  that 
this  model  has  been  successfully  employed  by  Gross  (31)  and  by 
Pitaevski  (32)  to  demonstrate  the  possibility  of  existence  of 
quantized  vortex  motion  In  a  Bose  fluid  and  to  Investigate  in¬ 
homogeneous  states  of  a  Bose  system.  We  employ  this  model  to 
investigate  the  nature  of  the  inhomogeneity  at  the  free  surface 
of  a  Bose  liquid  and  compute  thereby  the  energy  associated 
with  a  unit  area  of  the  surface.  This  leads  quite  naturally  to 
a  nonvanishing  value  of  the  intrinsic  surface  tension  ck(o), 
whence  follows  the  vibrational  part  of  the  surface  tension, 
a  (o),  so  that  finally  we  have  an  estimate  for  the  observed 
surface  tension  at  T=0°K.  The  resulting  theoretical  estimate 
of  0.28  erg  cm  compare's  favourably  with  the  experimental 
estimate  of  0.37  erg  cm-2(33^3). 
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3.2  THE  IMPERFECT  GAS  MODEL 

As  a  model,  we  consider  the  case  of  particles  interacting 
by  means  of  weak,  two-body  forces,  taking  the  self-consistent 
field  theory  as  the  first  approximation  (33).  We  shall  assume 
that  the  range  of  the  interparticle  force  is  much  greater 
than  the  average  interparticle  spacing,  so  that  the  rapidly 
changing  forces  on  a  particle  are  less  important  than  the 
smoothly  varying  average  force.  Under  these  assumptions  one 
can  start  the  analysis  with  the  Hartree  approximation. 

We  assume  the  ground  state  wave  function  to  be  a  product 
of  single-particle  wavefunctions . 

N 

y  (x  ,  ...,  x  •  t)  =  II  g(x.  ,t) 

—  ^  i  =  l 

1 1  g(x,t ) | 2d3x  =  1  ,  (3-1) 


where  g(x^,t)  is  a  normalized  single-particle  state.  We  also 
assume  that  g(x^,t)  is  separable  such  that 

g ( x  ,  t )  =  g(x)  exp 


The  single-particle  wave  function  g(x,t)  is  thus  an 
approximate  stationary  state  of  the  hamiltonian 


N  p  .  2  n 

H  =  E  ±  £  V(x.  -  x.)  . 

i  =  l  2m  2  i*J 


(3-2) 


If  we  renormalize  the  wavefunction  g(x)  and  let  iJj(x)  - 
/!T"g(x),  we  then  have  /l'i'(x)  |2d3x  =  N  and  |ip(x)|2  will  represent 
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the  density  of  the  fluid.  In  terms  of  the  time-dependent 

self-consistent  field  equation  becomes 


V2ip  +  ip  ( x ) 


V(x-y_)  d 3y . 


(3-3) 


The  energy  of  a  stationary  state  is 

r 


E  = 


2m 


Vip ( x )  |  2  d3x  +  4 


^(x)  I  V(x-y )  |  ip(y )  |  2d3xd3y  .  (3.4) 


The  complex  wave  function  ip(x,t)  thus  governs  the  behavior  of 
the  superfluid  condensate. 

The  analysis  can  be  much  simplified  if  we  replace 
V(x-y)  by  the  approximate  pseudo-potential  Vc6(x-y),  where  (34) 

v©  (  m  ;  • 

This  does  not  alter  the  general  nature  of  the  results  for 
phenomena  whose  characteristic  length  -  the  thermal  wave- 
length,  X  =  h/2TrmkT)  -  is  larger  than  the  range  of  the  force. 
The  simpler  equation  is  then 

V2  +  |iH2j  ip  (3.5) 

where  we  see  that  y ,  the  chemical  potential,  is  just  the  energy 
increase  in  the  system  when  one  more  particle  is  added  to  it. 

Depending  upon  the  boundary  conditions  used,  solutions  of 
the  equations  of  the  type  (3*5)  can  often  be  found  in  terms  of 
elliptic  functions.  Abrikosov  (35)  has  discussed  approximate 
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solutions  of  this  equation  in  relation  to  type  II  supercon¬ 
ductors  .  In  the  case  of  the  noninteracting  Bose  gas  contained 
in  a  box  of  macroscopic  dimensions  with  rigid  walls  the 
solution  possesses  the  absurd  feature  of  exhibiting  an  enormous 
density  near  the  center  of  the  container.  Solutions  of  Eq.(3*5) 
have  the  virtue  of  giving  an  almost  uniform  density  distribution 
in  the  system,  except  near  the  rigid  walls  of  the  box,  or  near 
the  vortex  cores.  Ginzburg  and  Pitaevski  (36)  have  used  a 
similar  equation  in  order  to  discuss  the  lambda  temperature 
shift  in  liquid  helium  films .  In  many  cases  it  is  useful  to 
apply  certain  integral  relations,  discussed  by  Schiff  (37), 
Derrick  (38)  and  by  Darewych  and  Schiff  (39),  which  can  con¬ 
siderably  simplify  integrals  of  the  type  (3-4). 
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3.3  SURFACE  NON-UNIFORMITY  AND  FREE  ENERGY 

We  shall  apply  Eq .  (3-5)  to  study  only  the  surface  part 
of  the  fluid,  for  it  is  in  this  part  alone  that  we  have  a 
variation  of  n(r)  from  the  bulk  value  to  a  vanishingly  small 
value,  over  a  region  characterized  by  some  sort  of  a  "healing 
length".  It  is,  however,  important  to  note  that  the  existence 
of  a  free  surface,  which  is  one  of  the  basic  characteristics  of 
a  liquid,  as  opposed  to  a  gas,  necessarily  requires  that  the 
atoms  in  the  surface  region  of  the  liquid  experience,  consequent 
upon  a  small  outward  displacement,  a  predominantly  attractive 
force  toward  the  interior  of  the  fluid;  consequently,  the 
effective  scattering  length  of  interaction  for  the  atoms  in 
this  so-called  surface  layer  must  be  negative  ( 4 0 )  ,  I.e.  a=-|a|. 

Now,  considering  a  free  surface  normal  to  the  x-axis,  so 
that  the  operator  V2  in  (3*5)  Is  just  d2/dx2,  and  denoting  the 
value  of  n(x)  at  the  point  of  inflexion  (see  Fig  3.1)  by  n  , 
we  obtain  for  the  parameter  y. 


U  = 


4  Trail 2 

m 


(3.6) 


The  differential  equation  (3-5)  can  then  be  written  as 


2 

^-4-  +  8tt  I  a  I  3  =  8tt  I  a  |  n  ijj 

dxz  1  IY  1  1  sr 


(3-7) 


Introducing  the  reduced  variables 


^ _(x)  =  r|;(x)/n  1//2;  £  =  x/Z  , 

O  o 


(3.8) 
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where  £ 


(  87m  I  a  I  ) 


-1/2 


the  differential  equation  becomes 


s 


(3.9) 


The  relevant  solution  of  this  equation  that  we  shall  consider 
is  the  one  that  vanishes  at  £>°°,  and  has  a  vanishing  slope  at, 
say  £,z o.  This  corresponds  to  a  density  varying  from  a  bulk 
value  n(o)  at  x=o  to  a  vanishingly  small  value  of  over  the 
distance  of  a  healing  length  £.  For  all  §<o,  i.e.  everywhere 
in  the  bulk  of  the  liquid,  we  shall  take  n(x)  =  n(o)  =  nQ . 

We  find  as  a  solution 


ib  =  /2~  sech  £ 
ro 


(3.10) 


and,  accordingly,  for  the  superfluid  density 


n;(x)  =  2ng  sech2  (j)  . 


(3.11) 


As  x-^o,  the  particle  density  approaches  the  value  2n  ,  which 

o 

should  be  equal  to  the  bulk  value  nQ .  Thus,  the  point  of 

inflexion  of  the  (ij>  ,£)  curve  (see  Fig  2.1)  corresponds  to  a 

particle  density  exactly  one-half  of  the  bulk  value,  which  is 

clearly  a  very  satisfactory  situation.  In  terms  of  the  bulk 

density  n  the  solution  is  thus 
J  o 


(3.12) 


This  density  variation  is  similar  to  that  found  by  Cahn 
and  Hilliard  (40)  in  the  case  of  a  classical  fluid.  Cahn  and 
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Hilliard  minimized  the  free  energy  in  a  region  of  non-uniform 
composition  near  the  critical  temperature  and  obtained  the 
density  variation 


n  (x ) 


r' 


tanfe  (j)  +  1 


(3.13) 


—  1/2 

where  £a  (T-T  )  .  The  authors  also  calculate  the  corres- 

3/2 

ponding  interfacial  free  energy  to  be  proportional  to  (T-T) 
Considering  again  Fig.  2-1,  the  dashed  line  in  the  graph 
indicates  the  spatial  region  which  would  be  occupied  by  the 
particles  now  constituting  the  non-uniform  surface  layer  if 
they  were  spread  out  uniformly  with  a  density  equal  to  nQ  . 

Thus  K  =  (4ir  a  n  )  is  a  measure  of  the  extent  of  the 
spatial  inhomogeneity  in  the  free  surface;  in  magnitude,  it 
is  comparable  to  the  lateral  dimension  of  the  "core"  of  the 
quantized  vortex! 

The  structural  energy  per  unit  area  associated  with  the 
non-uniformity  in  the  free  surface  can  be  determined  by 
calculating  the  difference  between  the  energy  per  unit  area 
of  the  surface  layer  as  it  is  and  that  of  the  corresponding 
layer  if  the  distribution  were  uniform.  One  thereby  has 
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It  is  difficult  to  find  a  theoretical  estimate  for  the  scattering 
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length  a.  However,  Girardeau  (41)  has  discussed  the  Imperfect 
Bose  gas  under  the  influence  of  an  attractive  interaction  and 
has  arrived  at  an  estimate  for  a^-2A.  If  we  use  this  value, 
we  obtain  for  the  surface  tension  of  helium  II, 

0^(0)  =  0.09  erg  cm-2  (3-15) 

The  existence  of  the  intrinsic  surface  tension  would  lead 
to  the  presence  of  ripplons  on  the  free  surface  which,  in  turn, 
would  give  rise  to  a  zero-point  free  energy.  Atkins  (4)  has 

carried  out  this  calculation  earlier  with  the  result 

« 

v  , 

av(o)  =  fc  h  V)g(V)dV 

o 

/8cti(o)  Y 

7 it/4  \  p6'  /  ,  (3.16) 

which,  on  substitution  from  (3*14)  becomes 

Sir  ^  no1l4|a|1V 

a  ( o )  =  7  /3~  m  .  (3-17) 

v 

Substituting  the  relevant  values  In  this  expression  we  obtain 

a  (o)  =  0.19  erg  cm-2  .  (3-18) 

Combining  ai(o)  and  av(o),  we  obtain  for  the  observed  surface 

tension  0(0)  of  He4  an  estimate  of  0.28  erg  cm  ;  the  corres- 

-2 

ponding  experimental  estimate  is  0.37  erg  cm 

Of  course,  on  the  basis  of  a  gaseous  model,  though  interacting, 
we  cannot  expect  to  understand  the  behaviour  of  a  liquid 
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quantitatively,  especially  when  the  numerical  value  used  for 

the  scattering  length  is  not  very  reliable.  In  fact.  If  we 

0 

took  a  value  for  the  scattering  length  a  =  -4A,  the  agreement 
with  experiment  would  be  excellent.  Qualitatively,  In  any 
case,  the  theory  developed  here  indeed  appears  to  be  a  step 
in  the  right  direction. 

Finally,  we  remark  that,  since  the  temperature  depen¬ 
dence  of  the  surface  tension  of  liquid  helium  is  very  well 
explained  up  to  approximately  1.5°K  by  the  surface  energy 
contribution  of  the  surface  waves,  the  intrinsic  surface 
tension  ck  may  not  be  expected  to  be  strongly  temperature- 
dependent  near  T=0°K.  One  also  expects  this  from  consideration 
of  the  third  law  of  thermodynamics .  Since  the  entropy 
S  =  -  (grj-r)^,  we  can  write  the  surface  part  of  the  entropy  as 


3  a 
3T , 


where  s  is  the  surface  area  of  the  system  and  a  its  surface 
tension.  Thus,  since  the  entropy  Is  zero  at  T=0°K,  clearly 

^  =  0  at  T=0° K . 
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3.4  SURFACE  TENSION  OF  HELIUM3 

3  4 

Since  helium  is  a  Fermi  liquid,  while  helium  Is  a  Bose 

liquid,  the  difference  in  statistics  would  seem  to  imply ^that 
no  one  theory  could  satisfactorily  explain  certain  phenomena 
occurring  in  both  liquids.  However,  if  statistical  effects 
were  unimportant  for  certain  phenomena,  for  example,  the 
surface  tension  of  the  liquid,  then,  perhaps  one  theory 
could  give  adequate  results  for  both  liquids.  In  fact.  If 
statistical  effects  are  neglected,  the  theory  developed  above 
can,  quite  successfully,  be  applied  to  liquid  He3. 

We  will  show  that  the  surface  tension  of  liquid  He3  can 
be  understood  along  similar  lines  as  that  of  He4,  although  we 
cannot,  of  course,  claim  a  strong  theoretical  basis  for  this 
treatment.  We  shall  use  as  our  basic  equation  Eq.(3.5)>  but 
cannot  speak  in  terms  of  ground-state  wave  functions,  but  will 
use  Eq.(3*5)  as  differential  equation  for  the  density  distrib¬ 
ution  of  the  fluid  and  calculate  the  free-energy  associated 
with  the  density  non-uniformity  at  the  surface. 

Since  the  zero-point  energy  of  the  atoms  plays  an  important 
part  in  determining  the  surface  energy  of  a  liquid  at  low 
temperatures,  liquid  He3  would  be  expected  to  have  a  surface 
tension  much  smaller  than  that  of  liquid  He4.  In  fact,  the 
surface  tension  of  He 3  at  T=0°K  can  be  extrapolated  to  a  value 
of  0.153  erg  cm-2. 

Lovejoy(42)  has  found  that  the  assumption  concerning  the 
existence  of  surface  waves  also  explains  reasonably  well  the 
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temperature  dependence  of  the  surface  tension  of  He3  up  to 
approximately  1.5°K.  He  also  had  to  invoke  an  intrinsic 
surface  tension  in  order  tor  support  the  quantized  wave  motion. 
This  intrinsic  surface  tension  can  also  be  reasonably  well 
explained  by  the  theory  of  Section  3.2.  We  shall  adopt  a 
similar  approach  as  used  by  Amit  and  Gross  (43),  in  considering 
Eq.(3-5)  not  as  a  wave  equation,  but  as  an  equation  to 
determine  the  distribution  of  fluid  density  over  the  surface 
region.  Thus  we  rewrite  Eq .  (3.5) 


(3.19) 


where  y  is  a  parameter  which  can  again  be  determined  at  the 
point  of  inflixion  of  the  Curve  (3.1)  and  we  have 


Prom  this  equation,  we  obtain  as  a  solution  just  Eq.(3-ll) 


We  also  rewrite  the  equation  for  the  energy  due  to  the  surface 
non-uniformity  in  terms  of  the  density 


(3.20) 


Now  since  the  boiling  point,  critical  point,  and  density 
of  He 3  are  appreciably  less  than  the  corresponding  values  for 
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He4  (see  Table  3-*»l),  the  molecules  of  He 3  are  farther  apart 
than  the  molecules  of  He4.  Consequently,  the  attractive 
interaction,  while  still  present  in  the  surface  region.  Is 
much  smaller  in  He 3  than  In  He4.  The  scattering  length  then, 
while  still  negative,  should  be  considerably  smaller  in 

o 

magnitude.  If  we  take  a  value  for  a  -  -1A  we  obtain  a.(o) 

=  0.04  erg  cm  and  cr  (o)  =  0.11  erg  cm  .  We  can  thus 
obtain  good  agreement  with  experiment  up  to  approximately  1.5°K. 
Table  311  (44)  Some  Equilibrium  Properties  of  Liquid  Helium-3 

and  Helium-4 :  Normal  Boiling  Point  (N.B.P);  Critical 
Temperature,  Tc;  Critical  Density  pc;  Density  pQ  at  0°K. 


N.B.P . 

T 

c 

^c 

^0 

He3 

3 .19 

3.33 

0.0413 

0.082 

He4 

i — i 

C\J 

.=r 

5.21 

0.0675 

0.146 

Atkins  (4)  and  Lovejoy  (42)  have  suggested  that  the  reason 
why  the  surface  tension  of  He4  and  He 3  depart  from  theoretical 
predictions  at  approximately  1.5°K  is  that  the  simple  two- 
dimensional  Debye  model  does  not  take  Into  account  sufficiently 
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the  higher  frequency  surface  waves.  Secondly,  at  higher 
temperatures,  the  wavelengths  of  the  surface  waves  become 
comparable  to  the  Interatomic  distance  so  that  the  energy- 
momentum  relationship  for  ripplons  will  be  affected.  A 
third  reason  is,  of  course,  that  as  the  critical  temperature 
is  approached,  the  intrinsic  surface  tension  itself  begins  to 
decrease  rapidly,  to  disappear  at  T  .  It  is  probable  that 
the  temperature  dependence  of  the  intrinsic  surface  tension 
is  thus  responsible,  for  the  departure  of  the  experimental 
results  from  the  theoretical  predictions  around  1.5°K. 
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of  the  fluid  if  the  density  were  uniform  throughout. 


CHAPTER  4 


SURFACE  TENSION  AT  THE  He 3 -He 4  INTERFACE 
INTRODUCTION 

Since  1946,  when  Franck  (46)  suggested  that  a  phase 
separation  of  the  two  isotopes  He3  and  He4  might  be  expected 
at  low  temperatures  due  to  the  fact  that  He3  would  not  take 
part  In  the  superfluid  phenomena  characteristic  of  liquid 
helium  II,  much  research  has  been  done  on  mixtures  of  He 3 
and  He  4 . 

Phase  separation  was  discovered  experimentally  by 
Walters  and  Fairbank  (47).  As  the  temperature  of  the  He3-He4 
mixture  is  lowered,  a  critical  temperature  Is  reached  at 
which  phase  separation  begins  to  occur.  As  the  temperature 
approaches  zero,  the  upper  phase  becomes  very  nearly  pure, 
the  He4  fraction  decreasingly  exponentially  with  temperature. 
However ,' contrary  to  initial  expectations,  the  concentration 
of  He 3  in  the  lower  phase  does  not  go  to  zero  as  the  temper¬ 
ature  of  the  mixture  approaches  zero.  This  finite  solubility 
of  He 3  in  He4  was  suggested  in  1961  by  van  Leuwen  and  Cohen(48) 
on  the  basis  of  a  hard  sphere  model.  Measurements  by 
Edwards  et  al  (49)  of  the  phase  separation  curve  by  means  of  a 

o 

heat  capacity  method  showed  a  limiting  solubility  of  He  in 
He4  in  the  lower  phase  of  6J  at  T  =  0°K. 

Liquid  He3-He4  systems  also  exhibit  superfluidity. 

However,  down  to  very  low  temperatures.  He 3  does  not  partake 
in  superfluid  flow,  but  forms  part  of  the  normal  fluid.  In 
fact,  when  the  temperature  is  suddenly  raised  in  a  region  of 
the  liquid  mixture,  there  will  be  a  flow  of  normal  fluid  and 
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He  away  from  this  region  and  a  corresponding  flow  of  super- 
Huid  toward  it.  Consequently,  the  equilibrium  distribution 
of  He3  in  the  mixture  is  very  sensitive  to  small  heat  currents, 
a  fact  which  strongly  affected  initial  experimental  results 
on  He3-He4  mixtures. 

In  1949,  Abraham,  Weinstock,  and  Osborne  (50)  discovered 
that  the  A-temperature  depended  upon  the  concentration  of 
He 3  in  He4.  This  was  contrary  to  expectations  based  on  the 
classical  ideal  solutions,  which  would  have  a  A-point  Inde¬ 
pendent  of  concentration.  It  was  found,  however,  that  the 
A-temperature  decreased  almost  linearly  with  the  concentration 
of  He3  in  He4.  The  simplest  theoretical  derivation  of  the 
A-temperature  is  He3-He4  mixtures  was  advanced  by  Heer  and 
Daunt  (51)  in  1951.  They  treated  the  He4  component  as  an 
ideal  gas,  quite  independent  of  the  presence  of  He3,  and  the 
transition  temperature  of  this  ideal  gas  is 

T  .  h»  7  h  \# 

A  1 2 . 6l2V  /  (4.1) 

v  m  /  ? 


where  Vm  =  N^v^  +  N^v^,  v^  and  v^  being  the  mean  volumes  per 

atom  in  pure  liquid  He 3  and  He4  respectively.  The  change  In 

the  A-temperature  is  therefore  a  consequence  in  the  change  in 

V  /Ni.  and  hence 
m  h 


W 


T 


A 


o 


(4.2) 
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where  is  the  mole  fraction  of  He  3  in  the  mixture.  This 

theory  predicts  the  initial  slope  of  the  T^  versus  x^  curve 

to  be  -1.9a  whereas  experimentally  the  slope  (3T,/8x~)  of 

A  3 

the  curve  is  approximately  -1.5  (52). 

Another  property  of  liquid  He3-He4  mixtures  that  depends 
uppn  the  relative  concentration  of  the  two  components  is  the 
surface  tension.  The  surface  tension  of  the  mixture  can  be 
expressed  approximately  in  terms  of  the  surface  tensions  of 
the  pure  liquids,  i.e. 

o  =  X  +  ( 1-x )  +  ag  ,  (4.3) 

where  ae  is  a  small  positive  "excess."  surface  tension  (53). 

In  this  chapter,  however,  we  shall  be  concerned  with  a 
discussion  of  the  interfacial  surface  tension  of  a  liquid 
He3-He4  system  at  T  =  0°K,  i.e.,  after  phase  separation  has 

t 

occurred.  The  predicted  value  for  the  interfacial  tension 
depends  critically  upon,  the  model  one  adopts  for  the  surface 
tension  of  a  quantum  liquid.  In  Section  4.2  we  will  calculate 
the  value  of  the  interfacial  surface  tension  as  predicted  by 
the  "classical"  theory  of' surface  tension,  which  has  been 
applied  by  Amit  and  Gross  to  liquid  He1",  and  in  Section  4.3 
we  will  calculate  the  value  of  the  Interfacial  tension  as 
predicted  by  Atkins’  model. 
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4.2  CLASSICAL  INTERFACIAL  SURFACE  TENSION 

In  the  "classical"  theory  of  the  surface  tension  of  a 
liquid,  the  zero  point  energy  of  the  surface  modes  of  vibration 
plays  no  part.  This  is,  no  doubt,  correct  in  the  case  of 
classical  liquids  in  which  case  the  zero  point  energy  of  surface 
modes  of  vibration  is  only  an  infinitesimal  part  of  the  total 
value  of  the  surface  energy.  However,  Amit  and  Gross  (54)  have 
used  the  classical  theory  to  calculate  the  surface  tension  of 
liquid  He4  and  obtain  reasonable  experimental  argument.  The 
authors  assumed  a  trial  density  pattern  for  the  liquid 


P ( x )  -  4po  a2 


x 

( x  2  +  a  2 ) 


xs;a 


p(x)  =  Po 


x^a  (4.4) 


where  x  is  the  perpendicular  distance  from  the  wall  and  a  is  a 
parameter  to  be  determined  by  minimizing  the  surface  energy 
obtained.  The  energy  can  be  expressed  as  an  integral  by  using 


Eq .  ( 3 . H  )  . 


E/S 


2m 


d*  + 


dx 


+  V  p 

o 


dx  , 


(4.5) 


where  V  =  /v(r)  d3r.  Using  (4.4)  and  performing  the  integral 
one  arrives  at : 


E/S  =  mpQ 


0.73 


k2 

m2a 


0.16 


(4.6) 
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When  this  is  minimized  with  respect  to  a  one  obtains 

E/S  =  0.7  'h  c  p  0.4  erg  cm  2  ,  (4.7) 

which  agrees  reasonably  well  with  experiment. 

In  order  to  obtain  an  estimate  of  the  interfacial  surface 
tension  of  a  liquid  He3-He*  system,  based  on  the  classical 
model,  we  shall  apply  Antonoff's  rule  (55)  to  calculate  the 
interfacial  tension  0^(0)  at  T  =  0°K, 

0^4(0)  ~  ^4(0)  -  a^(o)  ,  (4.8) 

Cjj(o)  being  the  equilibrium  surface  tension  between  liquid  He4 
and  its  vapour  at  T  =  0°K  and  0^(0)  the  surface 

tension  between  liquid  He 3  and  its  vapour  at  T  =  0°K.  In  a 
number  of  experiments  on  a  phase-separated  phenol-water  system, 
Antonoff  et  al  (56)  found  that  the  measured  interfacial  tension 
agreed  to  within  experimental  error  with  the  prediction  from 
Antonoff's  rule.  Of  course,  we  must  first  correct  cr^(o)  for 
the  fact  that,  on  phase-separation  at  T  =  0°K,  the  lower  phase 
may  still  contain  a  6%  admixture  of  He3.  Hence,  using  Eq . 
(4.3)3  neglecting  the  small  excess  surface  tension  we 
obtain  thenresult 

0^4(0)  «  0.36  -  0.15  =  0.21  erg  cm-2.  (4.9) 

On  the  other  hand,  if  the  ’’quantum"  theory  of  the  surface 
tension  of  a  quantum  liquid  is  correct,  then  the  Interfacial 
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surface  tension  o^(o)  must  be  regarded  as  made  up  of  two 
contributions,  i . e . 

a34  ( ° )  =  +  a34v^0^  5  (4.10) 

where,  as  usual,  is  the  intrinsic  interfacial  tension 

and  c^^Co)  is  the  contribution  made  by  the  capillary  waves 
at  the  interface.  Cleaply  enough,  it  is  then  only  a 3 4 -^ ( 0 ) 
which  can  be  obtained  by. applying  the  Antonoff  rule;  and  again 
taking  into  account  the  6%  concentration  of  He 3  in  He1*  at 
T  =  0°K,  we  obtain 

a3i)1  (° )  “  am(o)  -  a3i(o) 

=  0.134  -  0.04  =  0.094  erg  cm-2.  (4.11) 

The  contribution  a 3 ^ v  must  be  calculated  separately  from 
the  zero  point  energy  of  the  capillary  waves  at  the  interface. 
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(4.3)  VIBRATIONAL  MODES  AT  THE  INTE&EACE 

Consider  a  liquid  of  density  pf  and  depth  hf  flowing  with 
constant  velocity  V’  over  a  liquid  of  density  p  and  depth  h 
which  flows  with  constant  velocity  V. 

We  take  the  x-axis  In  the  Interface  separating  the  fluids, 
with  y  being  measured  vertically  upward  from  the  interface  at 
equilibrium.  It  can  be  shown  (57)  that  the  complex  potential 
for  the  lower  fluid  is  just 

w  =  -  (V-c)z  -  -  t  — .  G  ; — t-  cos  m  (z  +  i  h  ),  (4.12) 

smh  m  h  5 

where  c  is  the  speed  of  propagation  of  the  surface  waves,, 
m  =  2tt/X  ,  and  a  is  a  constant.  The  streamline  ¥  =  o  then 
becomes  q  =  a  sin  m  x,  where  n  is  the  vertical  displacement 
of  the  interface  from  equilibrium.  The  complex  potential  for 
the  upper  fluid  can  be  obtained  by  writing,  V',-hif,  for  V,h. 

The  speed  in  the  lower  liquid  Is  then  given  by 


q 


2 


p-  dw  =  (v  )2 

dz  dz 


2ma(v-c) 2 

sinh  m  h 


cosh  m  (y+h)  sin  m  x. 


and  the  speed  at  the  interface  (y=o)  is  then 


qQ  =  (V-c) 


H  - 


2m  q  coth  m  h 


J. 


( 4 .13a) 


and,  for  the  upper  liquid. 


6 


q1  =  (V’-c)2  jl  t  2  m  n  coth  m  h'/. 


( it .  13b ) 
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Bernoulli’s  equation  for  the  two  liquids  just  gives 

p’  +  2  P '  ‘Iq  2  +  P’  g  h  --const.  (4.14a) 

1  o 

p  +  2pqo+pgr)  =  const-  (4.14b) 


However,  in  the  presence  of  an  interfacial  surface  tension  a 
we  have  the  relation 

pf  -  p  =  a/R  =  -  a  Ip1  =  am2n-  (4.15) 

Hence  subtracting  (4.14B)  from  (4.l4a)  and  using  (4.15)  we  have 

!/2  P’q^2  -  i  pq2  +  g  n(pf-p)  -  am2n  =  const.  (4.16) 

Now,  substituting  q^2  and  q2  from  (4.13)  and  setting  the 
coefficient  of  n  equal  to  zero  (since  the  r.h.s.  of  (4.16)  is 
constant)  we  obtain 

m  p  (V-c)2  coth  m  h  +  m  p’(V’-c)2  coth  m  h’ 

=  g(p-P  f )  +  am2 .  (4.17) 


Letting  h  =  -h '  =  00  for  our  situation,  we  obtain  the  result 


2  _  Ag_  p-p  ’  ,  2ttp 

2tt  p+p’  X  (  p  +  p  ’  ) 


(4.18) 


As  usual,  we  shall  consider  only  short  wavelength  vibration^  so 

that  the  effect  of  gravity  can  be  neglected.  The  cutoff  frequency 

y  can  be  determined  in  the  same  way  as  for  a  pure  liquid.  It 
c 

Is  essentially  the  dif f erence , in  the  number  of  particles  per 
unit  area  In  the  adjacent  layers  on  opposite  sides  of  the 
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interface  that  are  responsible  for  the  restoring  force  when  a 
small  displacement  of  the  surface  is  made  (assuming  that  the 
interparticle  interaction  is  essentially  the  same  for  both 
fluids).  Hence  the  cutoff  frequency  vc  of  these  modes  is  deter¬ 
mined  by  requiring  that  their  total  number  per  unit  area  of  the 
interface  be  equal  to  the  difference  between  the  number  of 
molecules  per  unit  area  of  the  high  density  phase  on  the  one  hand 
and  the  low  density  phase  on  the  other.  It  turns  out  that 
■fc  -  4.3  x  1010  sec  1  for  the  interfacial  modes  which  may  be 
compared  with  the  values  9-5  x  1010  sec-1  for  liquid  He 3  and 
1.5  x  101 'sec  ^  for  liquid  He4.  Carrying  out  the  relevant  cal¬ 
culation  for  the  zero-point  energy  of  the  interfacial  waves.  In 
the  same  way  as  was  done  in  Section  (3*  ),  we  obtain 

a 3 4V ( o )  ~  0.007  erg  cm-2  ,  (4.19) 

which  is  extremely  small  In  comparison  with  the  Intrinsic 
surface  tension  (4.11).  Consequently,  the  theoretical  estimate 
for  the  He3-He4  interfacial  tension  would  be 

o^^(o)  -  0.10  erg  cm-2  ,  (4.20) 

which  is  about  one-half  of  the  classical  estimate.  An  extra¬ 
polation  of  experimental  determinations  of  the  He3-He4  interfacial 
tension  should  therefore  be  able  to  decide  the  Issue  between  the 
two  estimates  and  provide  a  practical  justification  for  expressing 
the  surface  tension  of  a  quantum  '..liquid  as  a  sum  of  two  con¬ 


tributions  of  a  different  nature. 
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CHAPTER  5 

LAMBDA  TEMPERATURE  DEPRESSION  IN  HELIUM  FILMS 
5 . 1  INTRODUCTION 

Although  film  flow  had  been  observed  much  earlier,  the 

existence  of  a  helium  film  in  equilibrium  covering  a  vertical 
surface  in  contact  with  a  helium  II  bath  was  first  postulated 
by  Rollin  (6l)  in  1936,  to  explain  the  anomalous  character  of 
the  heat  flow  in  the  liquid.  The  existence  and  properties  of 
such  a  film  were  first  clearly  elucidated  in  a  series  of  ex¬ 
periments  performed  by  Daunt  and  Mendelssohn  (62)  in  1939- 
Some  of  the  more  important  of  their  results  were: 

1)  Helium  II  creeps  in  the  form  of  a  film  over  all  solid 
surfaces  in  contact  with  the  bath  level,  as  long  as  the 
temperature  of  these  surfaces  is  well  below  the  A  point. 

2)  The  transfer  of  the  film  is  nearly  frictionless  and 
the  rate  of  transfer  is  in  the  first  approximation, 
independent  of  the  pressure  head*  or  the  length  of  the  path 
to  be  traversed. 

3)  The  superfluid  flows  up  a  temperature  gradient.  Such 
a  film  flow  is  thus  an  example  of  ’’ideal"  superfluidity 
since  no  corresponding  normal  fluid  flow  occurs  in  this 
manner . 

4)  Saturated  helium  II  films  (in  equilibrium  with  the 
saturated  vapor)  achieve  a  thickness  of  approximately 
300  A  at  a  height  of  1  cm.  above  the  surface  of  the  bulk 
liquid . 

Unsaturated  films  -  films  in  equilibrium  with  a  gas  at  a 
pressure  p  less  than  the  saturated  pressure  p  -  have  also  been 
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studied.  However,  as  pointed  out  by  Rice  and  Widom  (63),  there 
is  a  possible  connection  between  saturated  and  unsaturated 
films.  If  we  imagine  a  saturated  film  at  a  height  H  above  the 
bath  in  equilibrium  with  the  saturated  vapour,  then  the  pressure 
at  large  distances  from  the  wall  is 


P 


PQ  exp 


r_mgH, 
^  RT  ; 


Hence  an  unsaturated  film  at  a  relative  pressure  p/po  can  be 
regarded  as  having  identical  properties  as  a  saturated  film  in 
equilibrium  with  the  saturated  vapour  at  a  height  H  above  the 
bulk  liquid  given  by 

H  =  In  (p/p  ) . 
mg  o 


Hence,  when  we  discuss  helium  II  films  in  this  chapter,  we  shall 
not  differentiate  between  saturated  and  unsaturated  films . 

Since  at  low  temperatures  the  helium  atom  is  chemically 
inert  the  formation  of  the  film  is  clearly  an  example  of 
physical  rather  than  chemical  adsorption.  The  force  necessary 
for  holding  the  film  to  the  wall  must  be  the  van  der  Waals 
attraction  between  the  helium  atoms  and  the  wall. 

Another  important  property  of  helium  II  films  is  that  the 
A-t emperature  in  such  films  may  be  considerably  depressed  below 
the  bulk  value,  and  the  magnitude  of  this  depression  depends 
upon  the  film  thickness. 

Prederikse  (6*0  in  19^9  noticed  that  in  helium  II  films, 
the  specific  heat  maximum  becomes  smeared  out  and  occurs  at 
lower  temperatures  than  the  bulk  specific  heat  maximum.  It 
also  appeared  at  this  time  that,  for  a  given  film  thickness, 
the  onset  temperature  for  superflow  was  well  below  that  of  the 
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specific  heat  maximum.  Later  experiments  (65)  showed  that  in 
very  fine  channels  of  diameter  approximately  70A,  the  two 
temperatures  were  the  same.  It  has  been  concluded  experimen¬ 
tally  by  Brewer  et  al  (66)  that  the  disparity  between  the 
specific  heat  maximum  and  the  onset  temperature  for  superflow 
in  thin  films  is  consistent  with  the  presence  of  excitations 
associated  with  the  free  surface  of  the  helium  II  film.  As 
these  excitations  could  not  be  present  in  filled  fine  pores, 
the  two  temperatures  would  be  expected  to  agree.  These 
authors  conclude  that  for  a  given  film  thickness,  superfluidity 
dies  out  at  a  well-defined  temperature  above  which,  however, 
the  excitations  contributing  to  the  specific  heat  still  remain 
until  the  ordered  phase  is  completely  destroyed.  The  experi¬ 
mental  curves  for  the  temperature  at  which  the  specific  heat 
maximum  occurs  can  be  fitted  near  the  bulk  X-point  as 

.-14 


AT= 


0.24  x  10' 


'K. 


(t-3.2  x  10-8)2 

The  temperature  for  the  onset  of  superfluidity  obeys  approxi¬ 
mately  the  empirical  relation 

,m  0.92  X  HT11*  o„ 

AT  (t-3.2  x  10-®) 2  K' 


Here,  the  first  layer  of 'helium  II  near  the  wall  Is  considered 
to  be  solid.  However,  there  is  at  present  no  agreement  on  this 
point  (6 7);  hence  we  shall  not  make. this  assumption  in  our  cal¬ 
culations  . 

In  this  chapter  we  shall  attempt  to  explain  theoretically 
the  depression  of  the  X  “-temperature  on  the  basis  of  the  existence 
of  excitations  at  the  free  surface  and  In  the  region  of  non- 
uniform  density  near  the  wall.  In  section  5-2  we  shall  outline 
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the  attempts  based  on  the  ideal  Bose-Einstein  gas  theory  to 
explain  the  T^.  depression.  In  Section  5*3  we  review  the 
Ginzburg-Pitaevski  theory  which  is  based  on  an  analogy  with 
the  corresponding  Landau-Ginzburg  theory  of  superconductivity. 
In  Section  5*4*  we  discuss  several  modifications  of  the 
Ginzburg-Pitaevski  theory,  and  in  the  subsequent  section  we 
present  our  own  theoretical  calculations. 
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5-2  FINITE  THREE-DIMENSIONAL  BOSE -EINSTEIN  SYSTEMS 

Since  the  A-transition  occurs  only  in  He^  and  not  in 

3 

He  ,  it  appears  that  statistical  correlations  must  play  a 

fundamental  role  in  accounting  for  the  ^-transition  and  for 

h 

many  other  superfluid  properties  of  He  .  There  are  many 

reasons,  however,  why  Bose-Einstein  condensation  cannot 

account  even  qualitatively  for  the  behaviour  of  He^  near  the 

A-transition.  Although  Bose-Einstein  condensation  of  the 

ideal  gas  takes  place  at  a  temperature  comparable. to  the  A- 

temperature,  the  thermodynamic  behaviour  of  the  Ideal  gas 

near  the  transition  temperature  is  very  different  from  that 

4 

of  He  .  For  example,  whereas  experimentally  one  observes  a 

logarithmic  singularity  in  the  specific  heat  at  the  A-point, 

the  Bose-Einstein  gas  predicts  a  maximum  in  the  form  of  a 

cusp.  Moreover,  whereas  the  ideal  gas  predicts  a  Cv  a  T3^2 

behaviour  near  T=0°K,  the  experimentally  observed  temperature 

dependence  of  the  specific  heat  is  just  that  due  to  phonons, 

namely  Cv  a  T3.  These  discrepancies  are  certainly  not 

surprising  since  the  theory  of  the  ideal  Bose-Einstein  gas 

does  not  take  Into  account  atomic  Interactions .  It  is  clear 

that  one  has  to  take  into  account  not  only  the  effect  of  the 

4 

short-range  repulsive  forces  in  He  ,  but  also  the  weak 

attractive  forces  which  are  responsible  for  the  liquefaction 
4 

of  He  at  low  temperatures . 

A  number  of  attempts  have  also  been  made  to  explain  the 
depression  of  the  A-temperature  in  thin  helium  films  on  the 
basis  of  the  ideal-gas  model.  Osborne  (68)  showed  in  19^9 
that  although  the  two-dimensional  Bose-Einstein  gas  does  not 
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undergo  true  condensation,  it  is  possible  to  define  an  accumul¬ 
ation  temperature  at  which  a  rapid  accumulation  of  particles 
into  the  ground  state  occurs.  This  accumulation  temperature 
is  non-zero  for  a  finite  two-dimensional  system,  but  tends  to 
zero  as  the  size  of  the  system  is  increased  without  limit,  as 
we  have  already  shown  in  Section  2.2. 

It  was  also  shown  by  Ziman  (69)  that  a  finite  three- 
dimensional  ideal  Bose-Einstein  system  behaves  In  many  ways 
like  the  two-dimensional  system.  One  can  similarly  define  an 
accumulation  temperature,  which  is  non-zero  when  the  dimensions 
of  the  system  are  of  comparable  length,  but  goes  to  zero  when 
any  two  of  the  dimensions  are  increased  without  limit.  The 
perfect  gas  theory  thus  predicts  that  the  A-temperature  should 
depend  on  the  size  and  the  shape  of  the  container  in  a 
definite  way:  for  a  given  thickness  of  the  helium  film,  the 
lambda  temperature  should  go  to  zero  as  the  two  lateral 
dimensions  are  made  much  larger  than  the  thickness  of  the  film. 
This  is,  of  course,  in  complete  disagreement  with  experimental 
observations . 

In  1953,  Ziman  (69)  suggested  an  ideal-gas  model,  with  the 

aim  of  reconciling  the  discrepancies  between  the  experimental 

results  and  the  predictions  of  the  ideal  gas .  Ziman  suggested 

that  the  model  to  be  adopted  for  liquid  helium  was  made  up  of 

a  very  large  number  of  finite,  ideal  Bose-Einstein  cells  whose 

4 

size  is  determined  by  some  characteristic  length  of  He  .  These 
cells  were  to  have  in  common  the  same  temperature  and  the  same 
chemical  potential  but  were  to  be  otherwise  completely  inde¬ 
pendent.  Then  for  those  films  that  have  a  thickness  less 
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than  the  size  of  the  characteristic  cells  or  assemblies, 
finite  size  effects  should  be  observed.  The  characteristic 
extent  of  these  assemblies  were  chosen  by  Ziman  in  such  a 
way  that  one  obtains  the  best  numerical  correspondence 
between  the  variations  with  film  thickness  of  the  predicted 
accumulation  temperature  and  the  actual  temperature. 

Ziman’s  model  predicted  a  satisfactory  value  for  the 
thickness  of  a  saturated  film  below  the  A-point  and  also 
suggested  a  possible  explanation  of  the  magnitude  of  the 
characteristic  velocity  of  flow  through  superleaks.  However, 
the  model  does  not  make  satisfactory  predictions  of  the 
thermodynamic  behaviour  about  liquid  helium  in  the  neighbourhood 
of  the  A-point,  neither  in  thin  films  nor  in  bulk  samples. 

Goble  and  Trainor  (70)  have  improved  upon  the  numerical 
results  obtained  by  Ziman.  By  means  of  extensive  computer 
calculations,  they  evaluated  the  thermodynamic  quantities  per¬ 
taining  to  Ziman’s  model  by  carrying  out  explicitly  the  summations 
involved  rather  than  replacing  them  by  integrations,  as  is 
customarily  done.  They  also  discussed  the  condensation 
process  in  finite  systems  and  showed  that  small  Bose-Eins tein 
assemblies  retain  important  features  characteristic  of  the 
condensation  process.  However,  they  noted  that  a  suggestion  by 
Khorana  and  Douglass  (71),  namely  that  the  procedure  used  by 
London  (72)  for  determining  the  condensation  point  in  bulk 
systems  was  to  be  used  for  finite  systems  also,  is  misleading. 
Goble  and  Trainor  point  out  that  there  is  no  unique  procedure 
for  defining  the  condensation  temperature  in  finite  systems . 
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In  fact,  the  authors  considered  six  different  criteria  for 
evaluating  the  transition  temperature  in  a  finite  assembly  of 
size  LxLxD,  where  D<<L  (corresponding  to  a  thin  film).  All 
the  six  criteria  can  be  used  in  the  case  of  bulk  systems  also, 
in  which  case  one  obtains  the  same  transition  temperature. 
However,  for  thin  films,  different  criteria  lead  to  different 
results . 

For  completeness  and  for  possible: later  reference,  we 
reproduce  in  Table  5.1  the  numerical  results  obtained  by 
Goble  and  Trainor.  The  different  transition  temperatures  have 
the  following  meanings: 


1.  Tc  corresponds  to  the  maximum  value  of 


2.  TCp  corresponds  to  the  minimum  value  of 


3.  Tl  is  obtained  by  setting  Eo-y=0  in  the  expression  for 
the  total  number  of  particles, 

4.  T^  is  the  temperature  at  which  the  maximum  occupation 
of  the  first  excited  state  occurs, 

5.  T  is  obtained  by  extrapolating  the  No/N  vs.  T  curve, 

6.  T^  is  the  temperature  at  which  the  specific  heat 
attains  its  maximum. 

The  experimental  results  for  helium  films ,  which  are  given  in 
the  last  column  were  obtained  by  Long  and  Meyer  (73).  The 
temperature  values  have  been  renormalized  in  order  to  bring  the 
various  condensation  temperatures  for  bulk  samples  into  agreement 
with  the  observed  A-temperature  of  2.1K°K  in  liquid  helium. 
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Thickness 

D  in  A 

Renormalized  Temperatures  in 

°K 

Experimental 

Te 

T 

cp 

T 

2 

T" 

*Ti 

Te 

T1 

Value  of  T^ 

0 

0.65 

0.52 

- 

0.49 

0.54 

OO 

- 

10 

1.36 

1.17 

0.94 

1.13 

1.39 

3-24 

1.46 

25 

1.98 

1.79 

- 

1.80 

1.91 

- 

1.83 

60 

2.17 

2.09 

2*03 

2.01 

2.14 

2.38 

2.14 

750 

2.18 

2.18 

2.18 

2.18 

2.18 

2.18 

2.18 

Bulk 

2.18 

2.18  . 

2.18 

2.18 

2.18 

2.18 

2.18 

- 

The  transition  temperature  versus  film  thickness  results 
in  the  first  five  columns  are  in  qualitative  agreement  with  the 
experimental  results.  However,  whereas  experimentally  one 
observes  that  the  specific  heat  maxima  become  less  prominent 
and  shift  towards  lower  temperatures  as  the  film  thickness  is 
reduced,  Ziman's  model  predicts  a  shift  of  the  specific  heat 
maxima  towards  higher  temperatures  as  the  film  thickness  is 
decreased.  As  the  authors  point  out,  this  is  a  defect  shared 
by  all  ideal  gas  models.  It  appears  therefore  that  a  correct 
explanation  of  the  X-transition  in  thin  films  cannot  be  based 
on  the  ideal  gas  model  but  needs  to  take  Into  account  the  atomic 
interactions.  Moreover,  it  is  possible  that  finite  size  effects 
alone,  as  considered  by  Goble  and  Trainor,  may  not  be  sufficient 
to  explain  the  X— temperature  depression  but  that  surface 
effects,  as  well  as  effects  due  to  the  van  der  Waals  Interaction 
between  the  helium  atoms  and  the  wall  have  to  be  taken  into 


account . 
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5.3  GINZBURG-PITAEVSKI  THEORY 


In  order  to  explain  the  phenomena  of  superconductivity 
and  superfluidity,  a  number  of  macroscopic  theories  which  do 
not  begin  with  the  elemental  particles  of  the  system,  but 
with  some  gross,  overall  features  of  the  whole  system,  have 
been  developed.  One  such  approach  was  developed  by  Landau  (7*0 
in  1937.  The  central  idea  of  that  approach  was  to  expand  the 
Gibbs  free  energy  of  the  system  in  terms  of  a  certain  order 
parameter  \jj  about  the  transition  temperature.  In  order  for 
this  approach  to  be  valid,  one  considers  only  the  region  in 
which  the  order  parameter  is  small,  i.e.,  the  region  near  the 
transition  temperature:  | T— T  | < <T  . 

In  1950  Landau  and  Ginzburg  (75)  applied  the  theory  to 
the  superconducting  phase  transition  with  great  success. 

Eight  years  later,  Ginzburg  and  Pitaevski  applied  the  same 
theory  to  the  superfluid  transition,  though  not  with  the  same 
measure  of  success. 

In  this  section  we  review  the  Ginzburg-Pitaevski  theory 
and  its  application  to  the  depression  of  the  lambda  temperature 
in  helium  films.  In  Section  5-3  we  shall  consider  some  useful 
modifications  of  this'  theory. 

The  first  assumption  made  in  the  theory  is  that  the  order 
function  \jj(r)  corresponds  to  an  effective  wave  function  of  the 
superfluid  component.  Since  ip( r)  is  assumed  to  be  a  complex 
wave  function,  we  can  write 

\fj(r)  =  x(r)  (5-1) 

The  real  functions  x(r)  and  $(' r)  defined  by  (5.1)  correspond  to 
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the  superfluid  mass  density  p^(r)  and  the  superfluid  velocity 
v  (r)  by  the  relations 

S 

P$  ( r )  =  m  x2  (r )  > 

vs(r)  =  |  Vf(r) .  (5.2) 

Another  assumption  made  in  this  theory  concerns  the  functional 
form  of  the  Gibbs  free  energy  g (\p)  .  Considering  a  system  in 
which  the  normal  component  is  at  rest  and  is  at  a  temperature 
close  to  the  A-temperature ,  we  can  expand  g(ip)  in  powers  of 
the  parameter  ip(r),  which,  of  necessity,  is  very  small. 

The  functional  form  for  g(\|>)  that  Ginzburg  and  Pitaevski 
take  is  the  following: 

g(ip)  =  gx  2+6/2|<t>r-  |  ViJ>  |  2  .  (5.3) 

Here  the  term  g^  is  the  Gibbs  free  energy  for  the  normal 
component;  the  next  two  terms  are  the  free  energy  contributions 
of  the  superfluid  component,  and  the  last  term  represents  the 
kinetic  energy  contribution  from  the  superfluid  component. 

We  then  minimize  g(ip)  with  respect  to  ip  in  order  to  deter¬ 
mine  \p( r)  from  the  resulting  Euler  equation  which  we  can  write 
as 

“!m  v2^_ail;  +  3|Y|2^  =  0-  (5**0 

t 

This  is  the  Ginzburg-Pitaevski  equation  which  has  the  same  form 
as  Equation  (3-5)  although  the  coefficients  have  completely 
different  meanings . 

The  third  assumption  concerns  the  pressure  and  temperature 
dependence  of  the  parameters  ot(p,T)  and  6(P,T).  Since  the 
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transition  is  characterized  by  the  solution  ^=0  above  and 
^ ¥■  0  below,  the  parameter  a(P,T)  must  vanish  at  the  transition. 
Also,  for  g (ip)  to  be  a  minimum  at  the  transition  point,  the 
coefficient  3(P,T)  of  the  fourth  order  term  must  be  positive 
at  the  transition.  The  transition  line  is  therefore  charac¬ 
terized  by  the  equation  a(P,T)=0  (76). 

Considering  the  transition  for  a  given  value  of  the 
pressure,  Ginzburg  and  Pitaevski  assume  (after  Landau  (7*0) 


a(T)  =  a'  (T-Tk  )  =  -a1  AT. 


(5.5) 


The  coefficient  3(T)  is  simply  put  equal  to  the  constant  3(TO  . 

There  is  no  reason  why  this  expansion  of  the  parameters 
tt(T)  and  3(T)  should  be  valid.  In  fact,  we  shall  see  in 
Sec.  (5.^0  that  the  above  expansion  leads  to  an  incorrect  tem¬ 
perature  dependence  of  the  superfluid  density. 

The  Ginzburg-Pitaevski  theory  can  be  applied  to  the  problem 

A 

of  the  transition  temperature  depression  in  He  films.  One 
first  simplifies  Equations  (5.3)  and  (5.^0  by  introducing  the 
dimensionless  parameters 


£  =  Jl'/(2ma)'1/2 


(5.6) 


The  Gibbs  free  energy  then  takes  the  form 


(5.7) 


and  the  corresponding  Euler  equation  becomes 


(5.8) 


For  a  film  of  thickness  L  in  the  z-direction,  say,  and 
infinite  in  extent  in  the  x  and  y  directions,  Ginzburg  and 
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Pitaevskl  took  the  boundary  conditions  corresponding  to  an 
impenetrable  wall,  i.e., 

^0(£).=  0  for  £=0  and  Z=L/l,  _  (5.9) 

These  boundary  conditions  are,  strictly  speaking,  valid  only 
for  a  film  enclosed  on  both  sides  by  impenetrable  walls.  The 
solution  of  Equation  (5.8)  with  the  boundary  conditions  (5.9) 
can  be  expressed  in  terms  of  an  elliptic  integral  of  the  first 


kind  (77) 


r* 


f 


}l>. 


dt 


'o  v7t  4-2t  +  c  ,  (5.10) 

which  can  be  solved  explicitly  in  terms  of  Jacobian  elliptic 
functions.  One  obtains 

/  K 


*o(«)  = 


sn 


(sTTt 1  y  kj 


/TTji  v/77pi 

where  the  modulus  M  is  given  by 


(5.11) 


j  =  2/1+X  K C)c)  (5,12) 

K(./0  being  the  complete  elliptic  integral  of  the  first  kind. 
The  constant  C  is  expressed  in  terms  of  Jb  as 

4-k2 

c  =  (T+I2 ) 2 


A  plot  of  the  function  ^  has  the  form  of  a  dome,  which  is 
symmetric  with  respect  to  the  middle  of  the  film  and  the  height 
of  which  decreases  with  decreasing  film  thickness.  Also, 
there  is  a  film  thickness  below  which  no  solution  with 
exists,  i.e.  the  film  ceases  to  be  superfluid.  The  order 
parameter  ^Q(£)  can  be  seen  from  Eq .  (5 .11)  to  vanish  when  k  =  0, 
and  hence  we  obtain  from  Eq .  (5.12) 


£  =  2k(0)  =  7T  , 


(5.13) 
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which  defines  a  film  transition  temperature  T  that  differs 

A 

from  the  bulk  transition  temperature  by  the  amount 

-7T2  t)2 
AT^  =  2m  a' 


IL 2 )  .  (5-14) 

The  constants  a;  and  3  can  be  evaluated  from  specific  heat 
data.  When  these  values  are  substituted  into  Eq ,  (5.14)  one 
obtains 


2x10  14 

“T7 


(5-15) 


One  finds  that  this  result  differs  from  the  experimental  values 
by  about  a  factor  of  3- 
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5.4  MODIFICATIONS  OF  THE  GINZBURG-PITAEVSKI  THEORY 

The  GInzburg-PitaevskI  theory  possesses  a  number  of  obvious 
shortcomings.  First  of  all,  as  has  been  mentioned  in  the 
previous  section,  there  is  no  basic  reason  why  the  expansion 
a(T)  =  a'  (T-TO  should  be  valid.  In  fact,  this  expansion  can 
be  shown  to  lead  Immediately  to  an  incorrect  temperature  de¬ 
pendence  for  the  superfluid  density: 

ma  mo/ 

ps  =  m  |  T/;e  |  2  =3  =3  (T-T^ )  ,  (5-16) 

which  disagrees  with  the  experimentally  well-established  relation 
*  (T^-T)27"3  near  the  A-transition  (37^38). 

Secondly,  the  Ginzburg-Pitaevski  theory  predicts  a  second- 

order  phase  transition  in  a  bulk  system.  The  Gibbs  free 

energy  in  the  bulk  limit  can  be  written  from  Eq .  (5-3) >  using 

the  relation  a(T)  =  a7  (T-T.  )  , 

a'  2 

g  =  gl  "  W~  *  (5.17) 

This  clearly  represents  a  second-order  phase  transition 

in  contrast  to  the  logarithmic  behaviq'ur  of  Cp  near  the  lambda 
point . 

Thirdly,  the  Ginzburg-Pitaevski  theory  predicts  a  second- 
order  phase  transition  in  films  also  (80).  Measurements  of  Cp, 
however,  show  a  very  smooth  behaviour  at  the  superfluid 
transition  temperature  indicating  thereby  a  higher  order 
transition  (81,82). 

There  are  a  number  of  obvious  ways  of  "improving”  the 
Ginzburg-Pitaevski  theory  to  obtain  better  agreement  with  exper¬ 
iment.  One  can,  for  example,  easily  Incorporate  the  log¬ 
arithmic  behaviour  of  Cp  near  T^  into  the  theory  by  keeping  the 
same  functional  form  for  g (VO  but  changing  the  temperature 
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dependence  of  the  coefficients  a  and  $  (80,83).  According  to 
such  a  scheme  one  chooses  a  =  -a"AT8(T),  S(T)  =  8"ln(-AT),  so 

i)  2  =  -a "AT, 

ii)  g  =  g1  -  |  a " 2  8 " ( AT ) 2  ln(-AT) .  (5.18) 

Consequently  one  obtains  the  proper  logarithmic  behaviour  for 
C  near  T, .  However,  the  theory  not  only  predicts  a  linear 

p  A 

dependence  of  superfluid  mass  density  on  temperature,  but 
when  applied  to  the  helium  film  shows  no  shift  at  all  in  the 
transition  temperature. 

Alternatively,  one  can  choose  the  temperature  dependence 
of  a  and  8  so  as  to  obtain  the  correct  relationship  between  the 
superfluid  density  and  the  temperature.  This  was  done  by 
Mamaladze  (84)  in  1966.  Mamaladze  assumed  that 


o  =  +a'(Tx  -  T)4/3,  e  =  B'(TX  -  T)2/3,  (5.19) 

where  a  and  8  are  both  positive.  Substituting  these  parameters 
into  Eq .  (5.16),  we  readily  obtain 

ps  =  W"  (TA  “  T)2//3  *  (5.20) 

We  thus  obtain  the  proper  dependence  of  superfluid  density  on 

temperature.  However,  on  substituting  the  parameters  a  and  8 
into  the  expression  (5.3)  for  g (ip)  we  obtain  the  Gibbs  free 
energy 
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which  again  implies  a  second-order  phase  transition  at  the  A- 
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(Tx  -  T)2, 


temperature . 

Modification  of  the  temperature  dependence  of  the  parameters 


a  and  8  also  leads  to  a  new  formula  for  the  A-temperature  shift 
in  films.  Substituting  Eq .  (5*19)  into  Eq .  (5-13)  we  obtain 
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ATa  =  2'5  ^  (5.21) 

This  does  indeed  give  a  better  agreement  with  the  experimental 
results  near  the  lambda  point.  However,  as  has  been  pointed 
out  by  Wong  (80)  amd  Mamaladze  (84)  there  is  as  yet  no  sound 
theoretical  basis  for  the  definition  or  redefinition  of  the 
parameters . 

Amit  (85*86)  has  further  modified  Mamaladze' s  theory. 

He  noticed  that,  whereas  the  requirement  that  \p=o  give  a  stable 
minimum  for  g(ip)  atvT=T^ ,  is. 'fulfilled  by  the  original 
Ginzburg-Pitaevski  theory  since  3,  the  coefficient  of  the  |  ip  j  4 
term  was  positive,  this  requirement  is  obviously  not  fulfilled 
by  Mamaladze 's  theory  since  both  the  coefficients  vanish  at 
T=T^  .  To  insure  stability  therefore  one  has  to  add  a  term  of 
the  form  y|*f|6  where  y>o.  We  then  have 

ski  =  g-L  +  (A+  |ravs2)f<2  +  BifiU  Cips,  (5.22) 

where  A  =  -a  (T^-T) 42  3 ,  B  =  gOl.-T) 2/3 ,  as  before.  Using  this 

theory  Amit  finds  for  the  critical  velocity  vc  in  a  film 

Ve  =  (2a/m)  (T^-T)2/3,  ,(5,23) 

and  for  the  superfluid  mass  density 

p(tv  =  0)  =  (T  -t)2/!  (5.24) 

s  A  - 

Wong  (80)  has  proposed  the  following  modification  of  the 

Ginzburg-Pitaevski  theory:  he  considers  an  order  function  ip , 

to  be  defined  later,  and  constructs  the  functional  g(ip)  as 

i  2 

g(iio  =  g;L  +  Aki2  +  Bkl3  +  ckP  +  ] vijj  1 2 .  (5.25) 

We  now  have  two  equations  to  define  the  lambda  point 
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By  varying  g(ijj)  with  respect  to  p,  we  now  obtain  a  sdmewhat 
different  Euler  equation 

-Jl2 

2m  V2ip  +  A \p  +  .  \jj 2  +  2 C\p3=  0  (5-27) 

From  microscopic  arguments,  Wong  is  able  to  show  that  the  order 

parameter  here  has  .a  different  meaning,  namely  that 

|  Tj>i  2  a  n  3/2  (5.28) 

s 

Consequently,  if  one  assumes  for  the  parameters  A  and  B  the 
temperature  dependence 
A(P,T)  =  At  At  , 

B(P,T)  =  BT  AT 

C (P ,T)  >  0  ,  (5.29) 

one  obtains  n  a  (-AT)2'3.  One  thus  obtains  the  proper  depen- 

s 

% 

dence  of  superfluid  density  on  temperature. 

Again,  one  obtains  a  second-order  phase  transition  for  a 
bulk  system.  However,  for  a  thin  film,  or  for  porous  media, 
Wong  obtains  a  third-order  phase  transition,  which  is  indeed  an 
improvement  over  the  original  Ginzburg-Pitaevskl  theory. 

Wong  has  also  applied  his  modified  theory  to  the  helium  II 
film  by  assuming  a  suitable  expression  for  the  order  function  ^ 
He  then  obtains 

(av  =  -  p w2°k  5  (5-30) 

which  is  in  very  good  agreement  with  experimental  results  near 
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5-5  DEPRESSION  IN  HELIUM  II  FILMS 

In  this  section  we  shall  present  results  of  our  calculations 
of  the  X— temperature  depression  in  films  using  four  different 
criteria.  This  will  tend  to  show  that  it  is  impossible  to  des¬ 
cribe  the  observed  depression  in  X-temperature  without  con¬ 
sidering  in  some  way  the  surface  effects.  The  calculations, 
though  their  merit  is  mainly  of  a  qualitative  rather  than 
quant ita-tioe  nature ,  will  show  that  these  surface  effects  become 
increasingly  important  as  the  thickness  of  the  film  is  decreased. 
Moreover,  the  calculations  are  found  to  show  a  reasonable 
agreement  with  the  experimental  data  concerning  the  specific 
heat  maxima  rather  than  with  the  data  related  to  the  breakdown 
of  superflow  in  films. 

The  first  criterion  that  we  use  for  determining  the  transition 
temperature  is  simply  to  equate  the  inertial  densities  of  all 
the  excitations  to  the  average  total  density  of  the  film;  for 

o 

the  latter  we  may  assume  the  bulk  .density  0.146  gm/cm'"  at 

a 

T=2.2  K.  We  know  that  this  density  varies  by  at  most  1%  between 
the  transition  temperature  and  absolute  zero  (22). 

In  order  to  calculate  the  inertial  density  of  a  gas  of 
excitations  we  shall- use  Landau's  definition  (87)  of  inertial 
densities.  According  to  this  method,  we  give  to  the  gas  of 
excitations  considered  a  drift  velocity  *v  in  a  certain  chosen 
direction,  and  calculate  the  resulting  drift  momentum  J.  The 
inertial  density  p  Is  then  taken  to  be 


...The  drift  momentum  of  a  gas  of  excitations  can  be  written  as 
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J  =  E  <n£>  p  =  K3 


N ( p )  p  d3pd3x  , 


(5.32) 


where  N(p)  represents  the  appropriate  occupation  number  for 
the  excitation  considered.  The  inertial  densities  of  the 
elementary  excitations  are  well-known  but  as  an  example  we  have 
calculated  in  Appendix  B  the  inertial  density  of  the  quantized 
surface  modes . 


The  inertial  density  of  rotons  turns  out  to  be 
4ttpo2  ( 27rykT)  e~A/kT  ^  ^ 


nr 


3h3kT 


pkT  +  p 


(5.33) 


For  phonons  we  obtain 

l67T5k4T4 

Pph  ~ 45h3c 3  ,  (5.3^) 

which,  at  temperatures  above  0.5°K  is  much  smaller  than  the 
roton  component. 

In  the  case  of  ripplons  we  obtain  the  result 

(^o)^  (kT)5,/3  r(j)  E(^)  (|)»  (5-35) 

which  depends,  of  course,  on  the  thickness  t  of  the  film. 

When  we  set  the  sum  of  the  inertial  densities  of  the  above 
excitations  equal  to  the  average  film  density  we  obtain  the 
results  given  in  Column  1  of  Table  5.2.  The  A-temperatures 
obtained  for  all  film  thicknesses  are  much  higher  than  the 
experimental  values  obtained  by  Long  and  Meyer.  Clearly  we 
have  not  taken  into  account  to  a  full  extent  all  the  surface 
factors  present  though  the  trend  seems  to  be  essentially  right. 

The  next  approach  is  somewhat  related  to  that  of  Ginzburg 
and  Pitaevski.  Since  the  free  energy  per  unit  volume  of 
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helium  II  is  a  continuous  function  of  temperature  even  at  the 
transition  point,  the  A-temperature  corresponds  to  a  definite 
value,  of  the  free  energy  per  unit  volume,  i.e.,  the  free 
energy  of  helium  II  must  reach  a  definite  value  before  the 
transition  to  the  phase  helium  I  occurs.  We  shall  assume  that 
this  value  of  the  free  energy  per  unit  volume  must  also  be 
attained  in  the  case  of  a  helium  II  film  in, order  that  one 
arrives  at  the  transition.  In  order  to  define  a  transition 
temperature  in  the  film  we  may,  therefore,  equate  the  free 
energy  per  unit  volume  in  the  bulk  system  at  T=2.l8°K  to  the 
free  energy  per  unit  volume  of  the  helium  II  film  at  T«T- . 

Since  surface  effects  in  the  film  contribute  considerably  to 
the  free  energy,  can  be  expected  to  be  considerably  modified 
in  relation  to  the  bulk  value  of  2.l8°K. 

For  the  time  being  we  shall  consider  the  free  energy  of 
the  film  to  be  made  up  of  three  parts,  namely  due  to  rotons, 
phonons  and  ripplons.  Needless  to  day.  It  Is  only  the  phonon 
and  roton  contributions  that  are  present  in  the  free  energy  of 
the  bulk  system. 

For  the  free  energy  of  phonons,  we  substitute  the  relevant 
energy-momentum  relation  E=pc  into  the  main  term  of  Eq.  (2.15), 
and  obtain 


-  8^£(4)_C mi  .  (5.36) 

In  the  case  of  rotons  also  we  have  to  calculate  the 


It 


volume  contribution  to  the  free  energy  in  a  similar  way  as  was 
done  for  the  surface  correction  In  2.3,  with  the  result  that 

Vi  kT  +  p02J  .  (5-37) 


(Dr  =  "  FT  ( 2itpk.T)  kT 


V2  ;  _  -A'/kT 
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Since  the  ripplon  surface  contribution  has  been  shown  to  be 
much  larger  than  the  surface  corrections  to  Eqs.  (5-35)  and 
(5.36),  we  can  neglect  these . 

The  ripplon  contribution  to  the  free  energy  of  the  film 
has  already  been  calculated  and  is  just  given  by  Eq .  (2.29), 

® a  =  -  ¥(^?)2/3  h(F)”  r(7/2)  p(7/3)  F  •  (5.38) 

Hence  we  have  to  solve  the  equation 


for  various  values  of  t.  The  results  are  given  in  Column  2  of 
Table  5.2,  and  are  in  somewhat  better  agreement  with  experiment 
than  those  in  Column  1.  It  is  clear,  however,  that  there  must 
be  other  factors  present  to  depress  the  A-temperature  even 
further.  This  bears  out  the  statement  of  Atkins  (89)  that  the 
lowering  of  the  A-point  in  thin  films  cannot  be  due  solely  to 
ripplons . 

However,  we  have  already  remarked  in  Chapter  3  also  that 
the  decrease  with  temperature  of  the  surface  tension  of  helium  II 
cannot  be  explained  solely  in  terms  of  the  ripplon  contribution, 
as  calculated  by  Atkins  (22).  This  may  be  due  to  two  reasons: 

1)  the  contribution  of  the  quantized  surface  waves  to  the 
free  energy  of  the  surface  may  not  have  been  taken  into 
account  properly,  or 

2)  there  are  additional  surface  excitations  present  which 
contribute  to  the  surface  free  energy. 
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Atkins  (22)  and  Lovejoy  (4l)  have  given  several  reasons  for 
believing  that  the  surface  wave  contribution  has  not 
sufficiently  been  taken  Into  account.  These  reasons  have  been 
reviewed  In  3.4.  Whichever  of  the  alternatives  is  true,  the 
total  surface  contribution  to  the  free  energy  can  be  taken 
into  account  by  taking  as  the  free  energy  due  to  surface 
excitations  the  experimentally  determined  result  of  Atkins  and 
Narahara  (21)  which  is  just 

/ F  )  _  -0.0081  T2-5 

\V/s  t  (5.40) 

If  we  substitute  this  into  equation  (5.38)  In  place  of  the 
ripplon  contribution  alone  (Eq.  5*40)  will  Include  the  ripplon 
component)  we  obtain  the  results  given  in  Column  3  of  Table  2.4. 
These  results  are  In  much  better  agreement  with  experiment. 

Next,  we  may  include  into  our  equation  a  free  energy  con¬ 
tribution  arising  from  excitations  near  the  wall.  Here  the  van 
der  Waals  attraction  between  the  helium  atoms  and  the  wall 
causes  a  nonuniformity  of  density,  as  shown  in  Appendix  II. 

We  have  presented  there  a  calculation  estimating  the  free 
energy  due  to  such  a  non-uniformity.  If  we  assume  that  the 
creation  of  excitations  with  increase  in  temperature  causes  the 
helium  film  to  be  less  loosely  bound  at  higher  temperatures, 
and  moreover  that  these  excitations  have  an  energy  spectrum 
similar  to  the  surface  excitations,  then  we  obtain  a  further 
contribution  to  the  free  energy  of  the  film  of  the  amount 


(5.41) 


0.0020  T2*5 


Taking  this  into  account  and  again  using  Eq .  5-38  as  our  basic 
equation,  we  obtain  the  results  shown  In  Column  4  of  Table  5*2. 
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Table  5.2 

Comparison  of  theoretical  estimates  of  the  ^.-depression  with 

A 

experimental  data. 


Thickness 

4  Itf  A 

Renormalized 

T  T 

1  12 

Temperatures 

T 

3 

in  °K 

Experimental 

Value  of  TA 

0 

0 

— 

0 

0 

- 

10 

2.05 

1.81 

1.66 

1.56 

1.46 

25 

2 . 11 

2.02 

1.97 

1.89 

1.83 

60 

2,16 

2.12 

2.09 

2.07 

2.14 

750 

2.18 

2.18 

2.18 

2.18 

2.18 

Bulk 

2.18 

2.18 

2.18 

2.18 

2.18 

The  results  appear  to  confirm  our  procedure  of  taking  into 
account  the  surface  effects  in  helium  films.  However,  there  are 
several  points  that  must  be  noted  here: 

1)  For  thick  films  the  effects  due  to  the  non-uniform 
regions  at  the  free  surface  and  near  the  wall  have  been  rather 
overstated.  When  the  film  is  thick  enough  the  central  layers 
of  the  film  can  be  considered  in  a  sense  independent  from  the 
"surface"  layers  and  hence  superfluidity  might  occur  in  the. 
bulk  of  fyh e  film  at  temperatures  higher  than  the  ones  calculated 
;in  column  4. 

2)  In  the  thin  film  region  the  experimental  values  lie 
below  the  calculated  ones.  One  reason  for  this  is  that  we 

have  treated  the  roton-like  excitations  as  mutually  noninteracting. 
It  is  clear,  however,  that  the  various  excitations  interact 
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strongly  at  temperatures  above  1  K.  Secondly,  and  equally 
significantly,  we  have  ignored  the  finite-size  effects  upon 
the  energy  spectrum  of  the  excitations  (47,90).  The  true 
A-point  depression  can  probably  only  be  explained  satisfactorily 
if  both  surface  effects  and  finite  size  effects  are  taken  into 
account  together. 

3)  The  measurable  "transition  temperatures"  in  helium 
films  are  (i)  the  onset  temperatures  for  superflow  and  (ii) 
the  temperature  at  which  the  specific  heat  of  the  system  is 
maximum.  Hence,  an  important  improvement  in  the  above  dis¬ 
cussion  would  be  a  calculation  of  the  specific  heat  maximum 
from  a  better  knowledge  of  the  energy  spectrum  of  the 
excitations . 

4)  One  could  also  do  the  calculations  considering  the 
first  layer  from  the  wall  to  be  solid,  as  done  by  Brewer  et 
at  (66),  in  which  case  agreement  with  experiment  would  be 
better  for  thin  films.  However,  as  there  is  as  yet  little 
agreement  on  this  question  ('67)  there  is  not  much  point  in 
doing  such  a  calculation  at  present. 

Nevertheless  it  appears  to  have  been  demonstrated  faihly 
conclusively  that  the  surface  effects  are  of  extreme  importance 
in  films  in  understanding  the  physical  behaviour  of  helium 


films . 
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CHAPTER  6 


DEPRESSION  OF  TRANSITION  TEMPERATURE  IN  ROTATING  HELIUM  II 
6.1  INTRODUCTION 

In  view  of  the  depression  of  the  transition  temperature  in 
helium  II  films  and  in  porous  media,  it  is  imperative  that  one 
consider  the  nature  of  the  superfluid  phase  transition  in 
rotating  helium  II.  A  depression  in  the  lambda  temperature  of 
rotating  helium  II  can  certainly  be  predicted  theoretically. 
However,  since  the  size  of  the  rotating  cylindrical  container 
is  characteristically  much  larger  than  the  corresponding 
thickness  of  the  helium  II  film,  it  will  be  evident  that  the 
mechanism  responsible  for  the  depression  of  the  transition 
temperature  in  rotating  He  II  differs  in  nature  from  the 
corresponding  mechanism  in  He  II  films .  For  the  same  reason  _ 
the  depression  in  transition  temperature  will  be  much  smaller 
in  the  case  of  the  rotating  He  II.  In  fact,  experiments  to 
this  date  have  not  yet  been  able  to  measure  the  lambda  temp¬ 
erature  of  helium  in  rotation  to  a  sufficiently  great  degree 
of  accuracy  in  order  to  detect  such  a  small  effect .  Before 
going  into  more  detail,  however,  we  shall  review  some  of  the 
relevant  properties  of  helium  II  in  a  rotating  medium. 

As  is  well-known,  when  a  cylindrical  vessel  containing 
helium  II  is  set  into  uniform  rotation,  the  whole  of  the 
fluid  appears  to  imitate  solid-body  rotation.  This  was  con¬ 
trary  to  the  prediction  of  the  Landau  theory  (88)  according  to 
which  the  .absence  of  viscosity  in  the  superfluid  component, 
and  the  necessity  of  maintaining  potential  motion  (i.e. 

curl  v  =0),  the  superfluid  component  would  not  partake  in 
s 
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the  rotational  motion  of  the  vessel  and  hence  only  the  normal 
component  of  the  liquid  would  rotate.  The  first  experiment 
to  test  this  prediction  was  performed  by  Osborne  (89)  in  1950. 
The  results  showed  that  helium  II  rotated  as  a  whole;  so  it 
was  thought  that  the  interaction  of  the  superfluid  component 
with  the  walls  of  the  vessel  or  with  the  normal  component 
must  be  responsible- for  this  "breakdown”  of  superfluidity. 
However,  an  earlier  experiment  performed  by  Andronikashvili  (90) 
in  19^6,  concerned  with  axial  oscillations  of  a  pile  of 
closely  spaced  parallel  thin  discs  had  shown  that  superfluidity 
was  not  destroyed  by  rotation  but  that  only  the  normal  com¬ 
ponent  was  dragged  into  rotation.  A  later  experiment 
(Andronikashvili  and  Kaverkin  (91))in  1955  showed  that  the 
fountain  effect  continued  to  exist  in  rotating  helium  II  and 
had,  in  fact,  the  same  quantitative  features  as  in  helium  II 
at  rest .  As  a  result  of  these  experiments  it  became  evident 
that,  although  an  interaction  did  exist  between  the  superfluid 
and  normal  components  of  rotating  helium  II,  this  interaction 
was  not  related  to  the  breakdown  of  superfluidity.  In  accor¬ 
dance  with  Onsager’s  hypothesis  (92)  in  19^9?  this  interaction 
was  related  to  the  appearance  of  vortex  lines  in  the  super¬ 
fluid  component  of  helium  II.  The  circulation  of  these 
vortices  was  quantized  in  units  of  h/m.  This  quantization  of 
the  circulation  explains  the  existence  of  an  initial  critical 

velocity  co  below  which  the  superfluid  component  does  not  at 
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all  participate  in  the  rotation.  For  if  a  vortex  line  cannot 
exist  in  a  state  with  arbitrarily  small  intensity,  then 
vortex  formation  can  take  place  only  if  some  minimum  energy 
is  available  in  the  fluid  flow.  Superfluidity,  on  the  other 
hand,  is  destroyed  only  when  some  greater  critical  velocity, 

say,  is  reached . 

c2 

Experimental  as  well  as  theoretical  evidence  for  the 
existence  of  vortices  in  rotating  helium  II  was  quickly 
forthcoming.  In  1955.?  Feynman  used  quantum  mechanical  ideas 
to  explain  the  possibility  of  vortex  formation  in  liquid 
helium  and  showed  that  the  established  fact  of  superfluid 
rotation  is  adequately  explained  by  the  existence  of  quantized 
vortices,.  The  macroscopic  velocity  field  formed  by  these 
vortices  for  each  velocity  of- rotation  simulates  that  of  a 
classical  fluid.  This  in  turn  defines  at  all  velocities  of 
rotation  an  equilibrium  density  of  vortex  lines  In  the  rotating 
fluid.  Experimental  evidence  for  the  existence  of  quantized 
circulation  came  from  the  classic  experiments  of  Vinen  (93) 
who  utilized  the  existence  of  a  Magnus  force  between  a  vortex 
line  and  a  vibrating  wire  to  demonstrate  that  circulation  did 
appear  in  quanta  of  h/m. 

The  first  experiment  designed  to  study  a  possible  shift 
in  transition  temperature  of  rotating  liquid  helium  was 
performed  by  Zamtaradze  and  Tsakadze  (94)  in  1964.  They 
studied  the  nature  of  the  thermomechanical  effect  near  the 
transition  point  In  order  to  determine  the  effect  of  a  uniform 
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angular  velocity  to  on  the  phase  transition.  They  found  that 
the  dependence  of  the  height  of  superfluid  helium  on  its 
temperature  does  not  differ  significantly  from  that  of  the 
fluid  at  rest . 

Bablidze,  Tsakadze  and  /Cheishvili  (95)  performed  an 
experiment,  again  In  1964,  which  was  designed  to  observe  the 
discontinuity  in  the  heating  or  cooling  curve  at  the  trans¬ 
ition  temperature  of  a  phosphor-bronze  thermometer  placed  in 
rotating  helium  and  then  compared  it  with  the  corresponding 
situation  for  liquid  helium  at  rest.  This  break  In  the 
heating  or  cooling  curve  is  connected  with  a  jump  taking 
place  in  the  thermal  resistivity  when  helium  II  passes  through 
the  transition  into  helium  I  and  vice  versa.  In  this  experiment, 
where  the  maximum  attainable  angular  velocity  was  to  =  0.54  sec  , 

no  shift  of  the  transition  temperature  to  within  the  experi- 

-4 

mental  accuracy  of  5  x  10  °K  was  observed. 

Theoretical  estimates  of  the  lambda  temperature  depression 
in  rotating  helium  II  have  been  made  on  the  basis  of  the 
Ginzburg-Pitaevski  model  of  helium  II.  KIknadze,  Mamaladze 
and  Cheishvili  (98)  have  considered  solutions  of  the  equation 

4-  U  -  V  +  w  x  r}2  vF-a'i'  +  3]'H2'i'  =  0  ,  (6.1) 

which  represents  a  generalization  of  the  Ginzburg-Pitaevski 
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equation  with  v^y  =  a  x  r,  i.e.,  the  normal  fluid  is  in  uniform 
rotation.  This  equation  admits  solutions  describing  a  unit 
vortex  line  in  an  Infinite  liquid.  As  x  -*•  0,  the  function  ¥ 
is  of  the  form  cx .  For  x>>£,  Eq .  (6.1),  with  63  =  0,  has  the 
asymptotic  solution 


where  £  =  ii  /  /2m a .  This  solution  represents  a  vortex  line 

In  which  the  density  at  the  vortex  core  equals  zero  and  rises 

to  the  normal  fluid  density  at  a  distance  of  the  order  of  the 

lealing  length  £.  Kiknadze,  Mamaladze,  and  Cheishvili  have 

shown  that  Eq .  (6.1)  also  has  solutions  corresponding  to  the 

presence  of  a  vortex  line  with  the  expected  difference  that 

the  vortex  axis  moves  with  the  normal  component  of  the  fluid. 

In  fact,  analogous  to  the  two-dimensional  array  of  magnetic 

field  lines  in  Type  II  superconductors,  as  determined  by 

Abrikosov  (96)  in  1957,  Eq .  (6.1)  also  admits  of  stable 

arrays  of  triangular  and  square  lattices  of  vortex  lines. 

The  possibility  of  the  destruction  of  superfluidity  Is 

realized  in  the  case  of  an  infinite  fluid  moving  with  velocity 

60  =  6)  .  As  60  60  ,  the  wavefunction  of  He  II,  obtained  as 

c  2  c2 

a  solution  of  Eq .  (6.1),  tends  to  zero  smoothly.  This  leads 
to  an  approximate  shift  in  the  transition  temperature  as  given 
by  (loj 

AT  =  -(2.3  x  10  1 1 ) Cjl>  °K 


(6.2) 
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where  w  is  measured  in  rad/sec. 

Kiknadze,  Mamaladze  and  Cheishvili  (98)  have  also  deter¬ 
mined  approximately  the  values  of  the  critical  velocity  of 
rotation  at  which  the  superfluid  density  in  a  vessel  of 
finite  size  becomes  zero.  One  thereby  obtains 

ATa  -  -(1.1  x  10"14)  /  R2  °K,  (6.3) 

where  R  is  measured  in  cm.  This  relation  is  valid  only  for  a 

small  range  of  angular  velocities  just  above  co  .  As  w  is 

C1 

further  increased,  AT^  will  depend  on  both  cd  and  R  until  the 

critical  velocity  cj  for  the  breakdown  of  superfluidity  is 

c2 

reached . 

In  the  next  two  sections,  we  present  a  calculation  of  the 
depression  of  the  transition  temperature  in  rotating  liquid 
helium.  The  functional  dependence  of  the  T, -depression  on  the 
angular  velocity  as  well  as  the  radius  of  the  (cylindrical) 
container  will  be  derived.  Further,  we  will  obtain  as  special 
cases  of  our  more  general  result  two  relations  which  are 
practically  the  same  as  relations  (6.2)  and  (6.3). 
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6.2  INERTIAL  DENSITY  OF  A  ROTATING  GAS  OF  EXCITATIONS 

Our  approach  consists  In  calculating  the  inertial  density 
of  a  rotating  gas  of  excitations  (phonons  and  rotons)  as  a 
function  of  the  angular  velocity  to,  the  radius  R  of  the 
cylindrical  container,  and  the  temperature  T.  This  will 
naturally  represent  the  density  p  (T)  of  the  normal  component 
of  the  fluid.  The  temperature  at  which  this  density  becomes 
equal  to  the  total  density  p  of  the  liquid  helium,  i.e.,  when 
pg  =  0,  we  will  take  as  the  transition  temperature.  Thus  a 
relationship  can  be  obtained  between  the  transition  temperature 
T^  and  the  parameters  to  and  R. 

We  employ  as  before  Landau’s  definition  (85)  of  the 
inertial  density  of  elementary  excitations  for  calculating 
the  inertial  density  of  a  rotating  gas  of  phonons  and  rotons. 

To  do  this,  we  study  the  gas  of  excitations  from  a  frame  of 
reference  In  which  it  appears  to  have  a  mass-motion  characterized 
by  a  drift  velocity  v  and  calculate  the  accompanying  drift 
momentum  J.  The  Inertial  density  p  Is  then  given  by  the 
obvious  formula 


.P  = 


jI  - 

v  V 


(6.4) 


where  V  is  the  volume  of  the  gas. 

The  drift  momentum  of  the  gas  can  be  written  as 


r 


J  =  £<n>p  =  r- 

-  n  £  -  h 


8  8 

n(p;v,m)  pd  xd  p. 


(6.5) 
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where  n(p;v,w)  is  the  relevant  mean  occupation  number.  For 
rotons  we  have 


n  ( p  :  v ,  oa ) 


-1 


-  exp  j'-8(Ep  -  v.p  -  to.  1)  j  ,  (6.6) 

where  E  =  A  +  (p-p0)2/2y,  and  6  =  1/kT,  and  £  is  the 

ir 

angular  momentum  of  an  excitation.  Carrying  out  the  Taylor 
expansion  of  (6.6)  about  (v,oo)  =  (0,0),  we  can  write 


2  3 

n(£;v,03)  =  n(p)  {1-66E  +  (6E)2  -  -|y  (6E)3  +...} 

(6.7) 


where  6E  =  -(v.p  +  o3.£_)  .  The  even-order  terms  of  this 
expansion  do  not  contribute  to  the  drift  momentum  (6.5) 
because  in  their  case  the  integrand  becomes  an  odd  function 
of  the  momenta.  The  first  two  odd-order  terms  yield  the 
following  result: 


J  =  N 


P  v 
kT 


1_ 

30 


P 


3 


v 


CkTT3 


{- 


|  02r2) 


J  . 


(6.8) 


where  R  is  the  radius  of  the  cylindrical  enclosure  and 
the  mean  number  of  rotons  in  the  gas  at  rest.  For  the 
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inertial  density  of  the  roton  gas,  we  thus  obtain 


Pr  (v,w)  =  pr  (0,0)  j\  +  yq 


P  2 
■  o 

kT 


P  R 
■  o 

20  skT 


(n£)~  v2  +  h  (^-)2u.2  +. 


.(6.9) 


It  may  be  pointed  out  here  that  the  first  two  terms  in  this 
expansion  are  just  the  first  two  terms  of  the  general  expan¬ 
sion  for  Pr(v)  as  given  by  Khalatnikov  (99): 


Pr  (v) 


Pr(0) 


(6.10) 


In  the  case  of  phonons ,  the  relevant  occupation  numbers 

are 


n(p ; v,w) 


(6.11) 


where  E  =  pc,  c  being  the  velocity  of  sound  in  the  liquid. 

J2. 

The  integration  (6.5)  is  once  again  straightforward,  though 
somewhat  laborious,  with  the  result 


(6.12) 


It  again  may  be  pointed  out  that  the  first  two  terms  of  this 
expansion  are  just  the  first  two  terms  of  the  general 
expression  for  Pph(v),  as  given  by  London  (100) 


P 


ph 


(v) 


pPh(0) 


(6.13) 
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In  the  limit  v+0 ,  we  get 

Pr(^)  =  Pr(0)  £l  +  w2  +  ...  J  ,  (6. 11!) 


and 


pPh(a,) 


pph(0) 


[■* 


R- 
c ' 


J 


(6.15) 


Since  we  shall  be  concerned  only  with  temperatures  in 

the  neighborhood  of  T=T^,  the  inertial  density  of  the 

phonons  can  be  shown  to  be  much  smaller  than  the  inertial 

density  of  the  rotons  and  may,  therefore,  be  neglected. 

However,  for  typical  values  of  the  parameters  involved,  viz. 

R  =  1cm,  and  w =30 rad/sec . ,  the  change  Ap^  in  inertial  density 

of  the  normal  fluid  can  become  quite  appreciable.  In  fact, 

for  ~  2.18°K,  we  have  a  relative  density  increase  as  high 

o 

as 


Ap 

r 

PrC0) 


1  x  10 


(6.16) 
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6.3  DEPRESSION  OF  TRANSITION  TEMPERATURE  IN  ROTATING  HELIUM  II 
Although  pycnometric  studies  by  Andronikashvili  and 
Tsakadze  (101)  had  apparently  shown  an  anomalous  increase  in 
the  density  of  helium  II  when  brought  Into  rotation,  of  a 
magnitude  comparable  to  that  given  by  Eq .  (6.16),  more  recent 
experiments  by  Esel’son  et  al  (102-105)  have  almost  con¬ 
clusively  proved  that,  when  sufficient  care  Is  taken  to 
ensure  temperature  stability  in  the  rotating  fluid,  one  does 
not  observe  any  density  changes  to  within  an  accuracy  of  1  in 
2  x  10^.  It  is  clear  then  that  Eq$,  (6.14)  and  (6.15) 
represent  an  increase  in  the  density  of  the  normal  fluid  at 
the  expense  of  the  superfluid  component.  Consequently,  the 
calculation  outlined  in  Section  6.2  is  valid  only  when  the 
superfluid  component  itself  is  in  rotation,  i.e.,  for 

03  ^  03  4  03  *  . 

C1  c2 

In  order  to  study  the  behaviour  of  the  transition 
temperature  T,  we  set  the  roton  inertial  density  (6.14)  equal 

A 

to  the  corresponding  quantity  for  non-rotating  helium  at  the 
transition  point  T^  .  We  then  obtain,  using  Eq .  (4.32), 


* 

s 


The  condition  0)^03  is  necessary 
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ituation  against  a  total  breakdown 
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-A/kT 
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-A 


T 
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o 
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o 


(6.17) 


Rearranging  (6.17)  in  order  to  solve  for  aT,(=T  -t1  L 

A  ' 


the  depression  in  the  lambda  temperature,  and  retaining  only 
the  first  two  terms  of  the  exponential,  we  obtain  for  the  T^ 
depression 


(6.18) 


AT 


X 


Substituting  for  the  relevant  numerical  factors  one  obtains 
the  final  result 


(6.19) 


where  co  and  R  are  measured  in  rad/sec.  and  cm  respectively. 

As  pointed  out  earlier,  this  formula  is  valid  for  the 
range  of  angular  velocities  in  which  the  superfluid  is 
dragged  into  rotation.  Blatt  and  Butler  (106)  have  shown  that 
the  first  internal  energy  jump  in  an  ideal  Bose  gas  contained 
in  a  cylindrical  vessel  of  radius  R  takes  place  at  an 
angular  velocity 


oo  =  4 . 45  ±  /mR2  . 
C1 


(6.20) 


This  is  comparable  in  value  to  the  critical  velocity  for  the 
creation  of  a  vortex  line  as  found  by  Arkhipov  (107)  and 
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Vinen  (93) : 

“c  *  I  (R2-a2)-‘  In  (5_),  (6.21) 

1  o 

where  aQ  is  the  core  radius  of  the  vortex.  Starting  from 
this  velocity  the  creation  of  vortex  lines  is  thermodynam¬ 
ically  favourable.  As  the  angular  velocity  increases,  the 
number  of  vortex  lines  increases  until  a  maximum  number  of 

vortex  lines  is  reached  at  an  angular  velocity  oo  ’  ^ 

C1 

(98). 

G1 

Hence  the  angular  velocities  for  which  the  formula  (6.19) 
for  the  ^-depression  is  valid  are  confined  to  a  narrow 
interval  of  angular  velocities 

4.45« 

w  >  mR*  .  (6.20)’ 

Substituting  this  value  of  co  into  equation  (6.19)  we  obtain 
for  the  dependence  of  AT^  on  R 

ATa  «  -(0.9  x  10~l4)  /R2  °K  ,  (6.22) 

which  is  in  excellent  agreement  with  the  corresponding  result 
of  Kiknadze  et  a  I  (12): 

ATa  -  -(1.1  x  10"14)  /R2  °K.  (6.22a) 


If,  on  the  other  hand,  we  decide  to  keep  the  radius  R  constant. 
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and  vary  the  angular  velocity  co,  then  we  can  obtain  the 
desired  result  by  eliminating  R  between  (6.19)  and  (6.20)’: 


ATX  ~  -(1.7  x  10  -11)  to  °K 


(6.23) 


which  is  again  in  good  agreement  with  the  corresponding 
result  of  Kiknadze  and  Mamaladze  (12): 


AT^  s  -(2.3  x  10  1X)  co  °K . 


(6.23a) 


One  must  indeed  remember  that  the  two  relations  (6,22)  and 
(6.23)  are  valid  only  for  constant  values  of  co  and  R  resp¬ 
ectively.  When  both  the  angular  velocity  and  the  size  of  the 
container  are  allowed  to  vary,  Eq .  (6.19)  must  clearly  be 
used  and  a  reliable  estimate  of  the  lambda  temperature 


depression  obtained. 
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CHAPTER  7 


THE  HELIUM  FILM  PROFILE 

7.1  INTRODUCTION 

A  large  number  of  experiments  have  been  performed  in  order 
to  determine  the  thickness  of  saturated  and  unsaturated  helium 
films  at  various  heights.  However,  due  to  the  difficulty  of 
guarding  against  small  amounts  of  radiation  influx  which 
greatly  influence  the  amount  of  evaporation  from  the  film, 
experimental  results  have  not  been  generally  in  agreement. 

Atkins  (108)  was  the  first  to  carry  out  a  systematic 
.measurement  of  the  thickness  of  a  saturated  helium  film  by  means 
of  a  dynamical  method  and  established  the  empirical  relation 

d  =  1.67  x  10-6/H1/7,  (7.1) 

where  both  d  and  H  are  measured  in  cm. 

In  a  series  of  experiments,  Burge  and  Jackson  (109)  and 
Jackson  and  Henshaw  (110)  have  found  the  empirical  result  to  be 

d  =  2  x  10“6/H1/3.  (7.2) 

However,  when  radiation  was  permitted  to  impinge  on  the  top  of 
the  film,  a  variation  approaching  that  of  Atkins  was  obtained. 

It  thus  appears  that  Atkins ?  result  was  not  reliable  due  to  the 
radiation  influx  occurring. 

An  experiment  by  Bowers  (111)  in  1953  showed  a  still 
different  variation  of  thickness  versus  height,  namely 
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In  a  similar  experiment  on  thin  unsaturated  films.  Bowers  found 
the  relation  (112) 

d  =  5  x  10~6/H1/3.  (7.4) 

Keeping  in  mind  the  possible  relation  between  saturated  and 
unsaturated  films,  as  pointed  out  by  Rice  and  Widom,  this 
result  is  reasonable  agreement  with  the  relation  obtained  by 
Jackson  et  al. 

The  earliest  theoretical  work  on  the  subject  was  done  by 
Schiff  (113)  and  Frenkel  (114),  who  took  into  account  the  nature 
of  the  van  der  Waals  force  of  attraction  between  the  wall  and 
the  adsorbed  helium  atoms  and  derived  the  relation 

d  =  (4  x  10"6)/H1/3  ,  (7-5) 

in  the  case  of  a  helium  film  on  a  glass  wall.  The  constant 
would,  of  course,  vary  with  the  material  of  the  Wall.  This 

3r 

V 

relation  is  in  good  agreement  with  the  results  of  Jackson  and 

coworkers  on  saturated  films  and  with  the  work  of  Bowers  on 

unsaturated  films  except  that  the  temperature  dependence  of 

1/2 

the  constant  k  in  the  relation  d  =  k/H  could  not  be 
determined  theoretically  whereas  the  experimental  value  varied 
considerably  with  temperature. 

A  different  theoretical  approach  was  adopted  by  Bijl,  de 
Boer,  and  Michels  (115)  who  considered  the  variation  of  the 
zero  point  energy  of  the  film  with  thickness.  They  visualized 
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the  eigenfunctions  of  lowest  energy  as  having  nodes  at  the  wall 
and  at  the  free  surface  of  the  film  and  having  wavelength 
X  =  2d,  thus  giving  a  zero  point  energy  per  atom 


(7.6) 


Minimizing  the  sum  of  the  zero  point  energy  and  the  gravitational 
energy  one  obtains  the  relation 


1  x  10  5 

H1/2 


(7.7) 


This  agrees  well  with  the  work  of  Bowers  on  thick  saturated 
films.  The  drawback  in  this  method  is  that,  with  the  eigen¬ 
functions  considered,  it  leads  to  a  concentration  of  helium 
atoms  toward  the  middle  of  the  film  which  is  physically  un¬ 
realistic.  However,  as  the  result  agrees  with  a  number  of 
experiments,  it  seems  that  the  correct  relation  between 
thickness  and  height  might  be  of  the  form 


d  =  kj/H1/3  +  k2/H1/2 


(7.8) 


In  this  chapter  we  shall  discuss  the  various  attempts  that 
have  been  made  to  derive  a  relation  of  the  form  (7.8).  In 
both  these  attempts  the  effect  of  the  finite  film  thickness  on 
the  zero  point  energy  of  the  helium  film  is  calculated,  and 
its  influence  on  the  d-H  relationship  studied.  A  comparison  of 
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the  results  obtained  by  different  approaches  makes  quite  an 
interesting  study. 

The  main  result  of  this  discussion,  however,  will  be  a 
vindication  of  the  contention  made  by  Pathria  (25),  namely 
that,  if  the  expression  for  the  density  of  single-particle 
states  is  written  as 

S(K)  =  S  v?  3  +  E(K)’ 

then  the  error  term  is  o(K2)  and  not  0(K2),  i.e.  the  error  term 
is  of  an  order  inferior  to  that  of  the  surface  term.  In  his 
paper,  Pathria  has  provided  computational  evidence  in  favor 
of  the  contention.  Here  we  shall  obtain  further  evidence  in 
favour  of  this  contention.  To  do  this,  we  will  present  a 
calculation  for  the  thickness  dependence  of  the  zero  point 
energy  of  the  film,  not  only  for  the  Neumann  boundary  conditions 
as  was  done  by  Atkins,  but  also  for  Dirichlet  and  periodic 
boundary  conditions.  The  results  will  show  that,  as  was 
contended  by  Pathria,  the  error  term  in  the  expression  for 
g(K)  is  definitely  Inferior  to  the  surface  term. 
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7 . 2  ZERO  POINT  ENERGY  OF  A  HELIUM  FILM 

After  Atkins,  we  assume  the  zero  point  energy  of  a  helium 
film  of  thickness  d  to  be  of  the  form 


(7.8) 


where  Z(d)  is  the  zero  point  energy  per  gram  of  a  film  of 
density  p  and  Z(°°)  is  the  zero  point  energy  per  gram  for  the 
bulk  liquid. 

Imagine  a  film  of  unit  area  at  a  height  H  above  the  bulk 
liquid  and  consider  an  increase  in  its  thickness  from  d  to 
d  +  t .  The  change  in  the  energy  of  the  system  is  then 


ptZ  (°°)  . 


In  equilibrium  AU  should  be  zero  and  hence 


(7-9) 


We  now  calculate  that  part  of  the  zero  point  energy  that 
arises  from  longitudinal  Debye  waves.  In  a  square  film  of  side 
L  cm.  and  thickness  d  cm.  a  typical  Debye  wave  of  frequency  v 


and  velocity  c  has  direction  cosines  1  ,  1  ,  and  1  giveh  by 
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nx>  ny 3  nz  being  integers.  Hence  the  lattice  points  corresponding 

to  waves  with  frequencies  not  exceeding  a  given  frequency  v 
will  lie  in  a  region  bounded  by  the  ellipsoid 


2  ,  2 
n  +  n 

*£ 


L‘ 


(7.10) 


The  length  L  of  a  side  can  be  made  very  large  permitting  large 

values  of  nx  and  n  3  but,  if  d  is  small,  the  maximum  possible 

value  for  n„  is  m  =  2d^c  where  v  is  the  cutoff  frequency.  We 
z  c  c 

choose  d  so  that  m  is  an  integer.  Then,  for  a  fixed  value  of  n 

Zj 


the  number  of  modes  with  frequency  between  v  and  v+dv  will  be 

2ttL^ 

g  ( n z ,  v  )  dv  =  — ~~  vdv  , 


but  with  a  minimum  value  v 


n  c 
z 

min  2d. 


The  cutoff  frequency  can  be  calculated  by  equating  the 
total  number  of  modes  of  vibration  to  N,  the  total  number  of 
particles  in  the  film. 
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for  Neumann  boundary  conditions,  and 

=  ~2  m  y2  -  ^  m(m+l)(2m+l) 
c 

for  Dirichlet  boundary  conditions,'.  Inserting  m  = 

00 

solving  for  \>  in  terms  of  v  ,  we  obtain 

c  c 


2d^ 

<+  and 


(7.11) 


where  \j  is  the  cutoff  frequency  for  the  bulk  liquid  and  the 
upper  sign  corresponds  to  Neumann  boundary  conditions. 

The  zero  point  energy  can  now  be  calculated. 


m 

L2  p  d  Z(d)  =  1 


v 

,  max 


1 

2 


h  v 


g(nz,v  )dv 


v  . 
min 


7rL2h 
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(m+l)v 


-  ^  m2 (m+1 ) 2 


for  Neumann  boundary  conditions. 
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and 


^  m2 (m+1 ) 2 


for  Dirichlet  boundary  conditions. 

O  H  < 

Substituting  for  m  =  ^c  ,  we  obtain  tor  the  zero  point 
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energy 
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and  using  relation  (7-11)  for  v  ,  we  have  finally 
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2(d)  =  Z(-)  [l  *  l  %)  +  \s  <f 


25  ,c 


ood 

c 
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...] 


(7.12) 


Substituting  this  expression  into  equation  (7-9)  we  obtain 


H  = 


a  _  ,  Z  (°°) 
gdr  32g 


) 


(7.13) 


irrespective  of  the  boundary  conditions  used.  Inserting  the 
relevant  numbers,  this  gives 


"  ■  ' 


(7.14) 


For  thin  films  this  agrees  with  the  experimental  findings  of 

Bowers  (112)  for  the  unsaturated  films  and  with  Jackson's 

results  (109 j  110)  for  the  saturated  films.  For  thick  films, 

however,  the  zero  point  energy  term  becomes  more  important  and 

1  /n 

hence  in  the  expression  d  =  k/H  ,  n  would  begin  to  approach  2. 
For  thick  films,  then,  this  result  agrees  fairly  well  with  the 
experimental  results  of  Bowers  (111)  on  saturated  films. 

In  the  case  of  periodic  boundary  conditions  the  wave- 

functions  will  have  the  form 
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where  k  =  (^-  ,  g-  ,  .  Again,  the  points  of  interest  are 

contained  in  the  ellipsoid 


Choosing  a  =  b  =  L,  this  corresponds  to  the  ellipsoid 


n  2  +  n  2  = 
x  y 


L2 


Again,  L  can  be  made  as  large  as  we  like,  but  the  maximum 

value  of  n  is  now  m  =  —  yr 
z  c  c 

The  number  of  modes  with  frequency  between  v  and  y+dy 

Q 

is  Aydy  with  the  minimum  value  of  y  being  ^  n  . 

Evaluating  the  cutoff  frequency  again,  we  obtain 
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Solving  for  v  in  terms  of  y°°,  we  have 
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We  now  calculate  the  zero  point  energy  as  before 


Z  ( d )  =  B  I 


n  =-m 
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max 
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m 


(7.17) 


Substituting  (7-16)  into  ( 7  - 1 7 (  we  obtain  the  final  result 


(7.18) 


This  is,  of  course  an  expected  result.  Finite  size  effects 
cannot  be  taken  into  account  properly  if  periodic  boundary 
conditions  are  adopted. 

We  can  use  the  foregoing  result  to  check  the  accuracy  of 
the  density  of  states  relation  (2.2).  If  we  do  the  calculation 
using  the  density  of  states  approach,  we  have 


g(?)dY  ~  td?  +  ... 


For  the  zero-point  energy  in  the  first  approximation  we  obtain 
in  a  straightforward  way 
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TT  Sh 
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Hence 

ZOd)  =  Z(»)  jl  ?  \  (^j)  +...J  ■  (7.21) 

This  result  Is  in  exact  agreement  with  Eq.(7-12);  hence 
E(K)  =  o(K2)  and  not  0(K2)  and  Pathria’s  contention  appears  to 
be  borne  out .  Still  further  evidence  for  the  correctness  of 
Pathria’s  contention  can  be  obtained  from  a  calculation  by 
Broy^ei- and  Pathria  (  US  ).  They  derive  an  expression  for  the 
free  energy  of  ideal  Bose-Eins tein  gas  particles  by  means  of  the 
density  of  states  formula  (2.2)  and  also  by  performing  the 
summation  over  states  directly  by  means  of  a  method  introduced 
by  Fowler  and  Jones  (118).  Complete  agreement  is  found  up  to 
the  surface  term.  Franchetti ( 117 )  and  Singh  and  Pathria(119) 
have  also  calculated  the  free  energy  of  ideal  gas  particles 
and  of  phonons  by  constructing  the  Mellin  transforms  of  the 
relevant  summands  and  have  obtained  complete  agreement  with 
the  results  obtained  by  using  the  density  of  states  formula. 
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In  the  case  of  periodic  boundary  conditions,  the  result 
obtained  by  using  the  density  of  state  formula  is  just  the 
first  term  of  Eq .  (7-21),  with  the  surface  term  being  zero. 
Hence  agreement  is  obtained  in  this  case  also. 

In  passing,  we  may  mention  that  Franchetti ( 117 )  has  also 
calculated  the  thickness  dependence  of  the  zero-point  energy 
of  a  film.  He  introduces  two  methods  for  calculating  the  zero 
point  energy  correction  in  films  and  his  result  differs  con¬ 
siderably  from  that  of  Atkins .  By  a  more  approximate  method 
Franchetti  obtains 


(7.19) 


Franchetti  does  not  obtain  a  term  proportional  to  1/d  because, 
when  summing  over  all  frequencies,  a  contour  is  chosen  so 
that  the  density  of  points  is  not  altered.  However,  Franchetti 
obtained  the  above  result  by  letting  the  total  number  of  modes 
of  vibration  equal  to  3N .  In  a  liquid,  however,  only  longi¬ 
tudinal  modes  of  vibration  are  present,  but  the  transverse 
modes  cannot  exist .  Hence  the  total  number  of  modes  of  vibration 
should  be  equal  to  just  the  total  number  of  particles  in  the 
liquid.  If  we  do  Franchetti' s  calculation,  by  calculating  the 
cutoff  frequency  from  the  condition  that  the  total  number  of 
modes  of  vibration  is  just  equal  to  N,  then  we  obtain  approximately 


(7.20) 


which  is  now  in  good  agreement  with  Eq .  (7-12). 
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CHAPTER  8 


CONCLUSION 

In  this  thesis  we  have,  in  the  first  place,  extended  the 
quantum  mechanical  theory  of  surface  tension  proposed  by 
Atkins,  which  considers  the  free  energy  contribution  of  quan¬ 
tized  surface  modes  of  vibration  as  an  essential  part  of  the 
surface  tension  of  a  quantum' liquid .  We  have  shown  moreover 
that  the  surface  correction  to  the  free  energy  of  bulk  liquid 
excitations  -  phonons  and  rotons  -  is  much  smaller  than  the 
free  energy  of  these  quantized  surface  modes,  or  ripplons . 

Our  main  contribution  to  this  theory  of  the  surface  tension 
of  liquid  helium  was,  however,  to  calculate  the  intrinsic 
surface  tension  of  liquid  helium  at  T  =  0°K.  This  contribution 
to  the  surface  tension  was  calculated  on  the  basis  of  the 
imperfect-gas  model  of  Gross  and  Pitaevski;  this  imperfect-gas 

model  is  able  to  give  a  description  of  the  density  non-unifor- 

surifc.e.e  ©P 

mity  at  thevliquid,  provided  one  assumes  an  attractive  inter¬ 
action  between  the  atoms  in  the  surface  region.  Consequently, 
one  calculates  the  free  energy  per  unit  area  associated  with 
this  non-uniformity  and  identifies  this  with  the  intrinsic 
surface  tension  of  the  liquid. 

A  test  for  the  validity  of  this  theory  of  the  surface 
tension  of  a  quantum  liquid  has  been  developed  in  Chapter  4. 

The  two  rival  theories  of  surface  tension  give  significantly 
different  predictions  for  the  interfacial  surface  tension  of  a 
phase-separated  He3-He4  system  at  T  =  0°K.  An  extrapolation  of 
experimental  measurements  on  the  interfacial  surface  tension  of 
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a  He3-He4  system  should  be  able  to  decide  the  issue  between  the 
two  estimates. 

In  Chapter  5  we  have  demonstrated  "qualitatively"  the 
necessity  of  taking  into  account  surface  phenomena  when  dis¬ 
cussing  the  lambda  temperature  depression  in  helium  films .  The 
criterion  that  we  have  used  to  define  a  lambda  temperature  was 
to  compare  the  free  energy  of  all  excitations  present  in  a 
helium  film  (including  surface  excitations)  to  the  free  energy 
of  the  excitations  in  bulk  liquid  helium.  Clearly,  since 
ripplons  contribute  to  the  free  energy  of  the  film,  a  lambda 
temperature  depression  will  be  observed,  which  depends  on  the 
film  thickness.  We  have  obtained  reasonable  agreement  with 
results  obtained  from  specific  heat  measurements,  although  the 
criterion  that  we  have  used  to  define  the  lambda  temperature  is 
not  one  that  can  be  directly  compared  to  experimental  results. 

In  Chapter  6  we  have  calculated  the  lambda  temperature 
depression  that  one  might  expect  in  rotating  helium  II.  Rather 
than  calculate  the  decrease  in  the  superfluid  density  as  a 
function  of  the  angular  velocity  and  the  size  of  the  rotating 
cylinder  as  was  done  approximately  by  Andronikashvili^,  _ 
Mamaladze  and  others,  we  have  calculated  the  increase  in  the 
normal  fluid  density  as  a  function  of  these  parameters;  this 
problem  can  be  solved  much  more  easily  and  surprisingly  enough, 
the  more  general  result  obtained  reduces  to  those  obtained  by 
Andronikashy  i'li  et  al  in  the  limit  in  which  the  lambda  temper¬ 
ature  depression  depends  on  either  the  angular  velocity  or  the 
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In  Chapter  7  we  have  presented  some  theoretical  justifi¬ 
cation  for  Pathria's  contention  that  the  error  term  in  the 
density  of  states  formula  (2.6)  is,  for  the  systems  that  we 
have  considered,  much  inferior  to  the  surface  term.  However, 
none  of  the  results  obtained  in  this  thesis  depend  critically 
upon  this  contention. 

The  main  criticism  to  be  made  of  the  above  thesis,  and 
of  much  of  the  previous  work  on  liquid  helium  is  that  the  results 
are  generally  based  on  macroscopic  theories  with  little  or  no 
basis  in  microscopic  theory.  The  main  reason  why  microscopic 
theories  have  not  been  as  successful  in  this  field  as  in  the 
field  of  superconductivity  is,  of  course,  that  helium  II  is  a 
liquid.  Consequently,  even  a  microscopic  theory  such  as  the 
hard-sphere  Bose  gas  cannot  be  expected  to  account  quantitatively 
or  even  qualitatively  for  many  phenomena  in  helium  II,  espec¬ 
ially  thqse  that  are  associated  with  the  free  surface,  or  with 
helium  II  films.  A  proper  microscopic  theory  must  take  into 
account  both  the  repulsive  and  attractive  nature  of  the  inter¬ 
particle  interaction  in  order  to  explain  without  ambiguity  the 
phenomena  occurring  in  liquid  helium  II. 
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APPENDIX  I 


FREE  ENERGY  OF  A  GAS  OF  NON-INTERACTING  BOSE  PARTICLES 
We  derive  here  expression  (2.28)  without  replacing  the 
summation  over  states  by  an  integration;  accordingly,  we  shall 
not  make  use  of  any  explicit  formula  for  g(K).  Let  us  consider 
a  gas  of  noninteracting  particles,  each  of  mass  m,  confined  to 
a  volume  V,  which  for  simplicity  is  taken  to  be  a  cube  of  side 
a.  The  energies  of  the  stationary  states  of  any  one  particle 
are  given  by 


e  =  (h2/8ma2)  (s2  +  t2  +  u2) 


(1-1) 


where  s,  t,  and  u  can  take  only  positive  integral  values  in  the 
case  of  the  Dirichlet  boundary  condition,  but  can  also  take  the 

value  zero  in  the  case  of  the  Neumann  boundary  condition. 

We  then  have  for  the  logarithm  of  the  grand  partition 
function 


(1-2) 


where  y=h2/( 8ma2kT)  and  z  is  the  fugacity  of  the  system.  The 
internal  energy  U  can  then  be  determined  by  the  relationship 


(1-3) 


where  $-=l/kT.  Carrying  out  the  expansion  (1-1)  in  powers  of  z 
and  the  differentiation  as  required  by  (1-2)  we  obtain 
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(1-4) 


The  sums  involved  here  can  be  evaluated  by  using  Poisson’s 
formula 
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and  the  one  resulting  from  (1-5)  by  a  differentiation  with 
respect  to  (yj),  namely. 
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From  (1-5)  and  (1-6)  it  follows  that,  in  the  limit  y-*0 , 
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Accordingly,  (1-4)  becomes 

u  =  I  (V/X3)kTg5/2(z)T  ^(A/A2)kTg2(z)  ,  (1-7) 

where  X=h/( 2irmkT)  'Z2  and  the  function 


n 


( z)  = 


00 

.  E, 
J=1 


•  -n  j 
J  zJ 


The  free  energy  of  the  system  can  be  calculated  by  means  of 
the  formula 


F  = 


U 

-  dT, 
T2 


(1-8) 


whence  we  finally  get 

F  =  (-V/Xs)kTg5/2  (z)  +  |  (A/X2)kTg2(z)  .  (1-9) 

For  T<T, ,  z-l;  then,  for  n>l,  g  (z)~C(n).  The  second  term  in 
A  n 

(1-9)  thereby  gives  precisely  the  expression  (2.28)  for  P  /A. 
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APPENDIX  II 


FREE  ENERGY -DUE  TO  NON-UNIFORMITY- AT  WALL 

In  order  to  calculate  the  free  energy  due  to  the  inter¬ 
action  between  the  helium  atoms  and  the  rigid  wall,  we  assume 
that,  due  to  the  short-range  van  der  Waals.  forces,  the  density 
distribution  of  the  film  as  a  function  of  the  distance  from 
the  wall  is  of  the  form  depicted  in  Figure  1. 


from  the  wall. 

r- 

At  present,  it  is-  impossible  to  take  into  account  the 
precise  nature  of  the  van  der  Waals  interaction  and  calculate 
the  density  nonuniformity  and  hence  the  free  energy  therefrom. 
Rather  than  take  into  account  this  interaction  directly,  let  us 
solve  a  problem  giving  a  similar  density  non-uniformity.  We 
shall  adopt  a  model  in  which  the  transition  region  I  is  a  region 
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of  density  variation  between  a  hypothetical  high  density  phase 
of  helium  II  and  a  low  density  phase  of  helium  II.  The  density 
of  the  liquid  in  region  I,  extending  in  the  negative  z  region, 
is  given  by  the  experimentally  determined  density  n  of  the 

3 

first  layer  of  helium  II  from  the  wall.  This  density  determines, 
or  is  determined  by,  the  nature  of  the  wall  and  the  van  der 
Waals  interaction;  hence  the  solution  of  this  hypothetical 
problem  determines  the  solution  of  the  actual  problem  of 
interest . 

We  describe  the  high  density  region  I  by  a  Gross- 
Pitaevski  equation 


(II-l) 


In  this  situation  we  let  y  represent  the  energy  change  that 
arises  when  one  particle  from  the  low  density  region  II  is 
transferred  to  the  high  density  region  I.  Hence  we  have 


=  -Malt2  (  n|2  =  -Ma|l2 

m  ' Y  i  1  m 


(n  -n  ) .  ( II-2 ) 

so 


Thus  the  square  of  the  wave  function  ip  represents  the  difference 

i 


in  density  between  region  I  and  region  II.  The  solution  of 
Eq .  (II-l)  can  be  obtained  analogous  to  the  solution  of  equation 
(3.9)  .  One  again  defined  dimensionless  parameters  and  obtains 


the  solution 
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where 


1 _ 

l  =  /4tt  |  a  |  (n  -n  J 
s  o 


At  any  point  in  the  transition  region  then  we  have 

n(x)  =  n  +  (n  -n  )  sech2(V^.  (II. 3) 

o  so 

If  we  compare  this  with  the  corresponding  solution  due  to 
Cahn  and  Hilliard  (3-22),  we  obtain  a  density  distribution 
which  is  more  or  less  similar.  These  authors  obtained  for  a 
transition  region  between  two  phases  of  a  classical  fluid 


n(x)  =  no  +  (ns-nQ)  tanh^  (n.D) 

Again,  as  before,  we  can  calculate  the  free  energy  due  to 
such  a  density  non-uniformity  and  obtain'  the  result 

F-F'  1  ^2(n  -n  ) 

o-  =  —  s  o 

A  '  J  m  X  (II. 5) 

In  order  to  obtain  a  numerical  estimate  for  this  energy  we 

must  again  estimate  the  value  of  the  scattering  length  a.  Since 

we  are  using  the  experimentally  determined  surface  tension 

0 

which  corresponds  to  a  scattering  length  a=-4A,  we  use  the  same 

value  here .  We  thus  obtain 

F-F  _p 

=  0.092  erg  cm  .  (II. 6) 

In  the  region  of  density  non-uniformity  near  the  wall, 
excitations  will  be  created  as  the  temperature  is  raised  above 
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0°K,  and  will  reduce  the  '’effective’'  free  energy.  If  we 
assume  that  the  temperature  variation  of  this  free  energy  is 
the  same  as  that  of  the  surface  free  energy,  we  obtain  for 
the  temperature  dependence  of  the  free  energy  due  to  the  non¬ 
uniformity  at  the  wall. 


F-F 

o 


2  R  -2 

erg  cm 


» 


A 


0.092  -  0.002  T 


(II. 7) 
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APPENDIX  III 


INERTIAL  DENSITY  OF  RIPPLONS 

In  order  to  calculate  the  inertial  density  of  a  two- 
dimensional  gas  of  ripplons  we  again  have  to  take  the  main 
term  in  the  expression  for  the  density  of  states  (2.11)  with 
n=2.  We  give  the  gas  a  drift  velocity  v  and  calculate  the 
drift  momentum  J.  The  drift  momentum  of  such  a  two-dimensional 
gas  of  excitations  is  thus 

J  =  J^n(e-v.p)  Pzd2p  .  (III-l) 

Applying  the  Taylor  series  expansion  (6.7)  for  n(e-v.]c), 
the  first  term  of  the  integral  for  the  Inertial  density  pn  of 
the  ripplons  becomes 

J  2  77  (  9n(  £  )  3j  /XTT 

pn  =  W  =  h*t  J  —  p  dp  •  (III-2) 

Using  the  ertergy-momentum  relation  e^  ?  =ffe)  V2  p3/2  for  ' 
ripplons,  and  substituting  into  Eq .  (III-2)  the  expression 

3n(  e )  _  2  /jip  1  5n(p) 

3e  3  \a  /  p 1/2  9p  , 


we  obtain 


5tt 

pn  3h^t 


3/2 
'  > 


(f2)  ^  J  n(p)  p 3/2 dp  . 

(k)^A  ’ 


(ni-3) 


we  substitute  for 


Letting  y  =  Ap?z,  where  A  =  ^ 
the  occupation  number  of  a  gas  of  ripplons  the  expression 
n(p)  =  (exp  Ap3/2-!)-1  and  obtain  the  result 


pn  = 


IOtt 

9h*t 


A-8/3 


y 


3/2 


ey-l 


dy  . 


This  expression  can  be  evaluated  in  terms  of  the  Riemann  Zeta 
function  to  yield 


n  v/3  • . 2  5/; 

pn  =  ^tW  r(5/3)  P ( 5/3 )  CkT) 


(I'll.  4) 
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